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(54) Multi-head nudear medldne camera for dual SPECT and PET Imaging 



(57) A nxitli-detector head nudear camera system 
automatically SMritchat)le (and optimized) to perform 
either SPECT imaging or PET imaging. The camera 
system employs, in one embodiment. muHi-detector 
configuration having dual head scintillation detectors 
but can be Implemented with more than tw» detector 
heads. The detectors contain switchable triggering cir- 
cuitry so that coincidence detection for PET imaging 
and non-coincidence detection for SPECT imaging is 
availat)le. Using a variat)le integration technique with 
programmable Integration interval, the event detection 
and acquisition circuitry off the camera system is switch' 
able to detect events of different energy distribution and 



count rate which are optimized for PET and SPECT 
Imaging. The system also includes dual integrators on 
each scintillation detector cfiannel for collecting more 
than one event per delector at a time for PET or SPECT 
mode. In PET or SPECT mode, the system also 
employs variat)le PMT cluster sizing having smaller 
cluster sizes for PET imaging an6 relatively larger dus- 
ter sizes for SPECT. In PET or SPECT mode, the sys- 
tem also employs variat)le centroid shape and zonal 
triggering. Utilizing the abwe programmable settings, 
the camera system can t>e automatically conf igured to 
operate in either SPECT or PET imaging modes. 
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Description 

BACKGROUND OF THE INVENTION 
5 (1) ReM Of the Invention 

The present invention relates to the field of nuclear medicine systems. Specifically, the present invention relates to 
signal processing systems for sdntillation detectors. 

10 (2) Prior Art 

Gamma cameras performing Single Photon Emission Computed Tomo^aphy (SPECT) have been utilized in 
nuclear medicine for some time. Anger proposed and developed such a system in the 1950s which has been modified 
and improved extensively with the introduction of high speed digital computer systems fbr ^ge acquisition as well as 

15 image reproduction. However. SPECT camera systems utilize a collimator tfiat is installed in front of the sdrrtillation 
crystal within a scintillation detector. The collimator is used to collimate tiie incoming gamma rays so ttiat only rays of a 
certain angle of incidence actually penetrate the crystal. Although SPECT imaging is extensively used in nuclear med- 
icine and provides ben^idal image quality, the coflimator introduces a source of image degradation in nuclear medicine 
inriages and tends to sornewhat reduce the resdution and quality of bnagesaoqu'^^ systene. 

20 Canoeras performing Posrlron Emission Tomography (PET) have been utilized in nuclear medicine as well witfi the 
introduction of relatively high speed detection electronics and computer systems for inr»ge acquisition and processing. 
These PET camera systems utilize a form of the scintillation detector that is used in SPECT systems, however, they do 
not utilize a collimator. In PET systems* the d^ection of two gamma rays in coincidence (in different scintillation detec- 
tors) is used to compute imaging information. A PET system employing two sdntillation detectors is descrbed n a 

25 paper presented by Gerd Muehllehner. M. R Buchin, and J. H. Dudek entitied "Performance Parameters of a Positron 
Imaging Camera," piMlshed in the IEEE Transactions on Nudear Sdence, Vdume NS-23, No. 1. on February 1976 
and also in a paper entitied "Performance Parameters of a Longitudinal Tomographic Positron Imaging System" by 
Paans, deGraaf, Welleweefd, Vaalburg and WoUring, in Nudear Instruments and Methods, Vdume 192, Nos 2.3, on 
February 1. 1982 pages 491*500. By utilizing higher energy gamma rays and diminafing the cdlimators, PETsystenrs 

30 offer greatiy improved image resdution and image quality over SPECT systems. Because the collimators are removed 
in PET systems, the detected count rate is higher in PET cameras over SPECT camera systems. Although both camera 
systems utilize different detection dectronics and other drcuitry, both PET and SPECT systems empicy sdntillation 
detectors. 

The detection hardware for SPECT and PET systems is different in temns of the manner in which the systems 

35 detect and record events and also different because PET systens operate at higher count rates over SPECT sys- 
tems. Further. SPECT'systems ennploy a different radionudide over PET systems and dated gamma rays at different 
energy levds over PET systems. For this reason, although SPECT and PET systems are versatile and useful within 
nudear medidne, in ttie prior art, different camera systems have been implemented and supplied for PET and SPECT 
imaging. Therefore, a facility desiring to perform SPECT and PET imaging is required to acquire two separate camera 

40 systems at a rdatively greater expense. 

tt would be advantageous, tiien, to provide a nudear camera system offering the ability to perform both SPECT and 
PET Imaging techniques within a single configurable system. Therefore, the expense of acquiring two separate systens 
can be advantageously avoided. The present invention offers such advantageous capabilrties. 

Accordingly, it is an object of the present invention to provide a single dual head nudear camera system that is 

45 swrtchaUe between SPECT imaging and PET imaging to p^form dtiier mode of operation. It is an ok)ject of the present 
invention to provide such a system employing high count rate event detection drcuitry for PET studies and also switch- 
able to offer detection of gamma rays of different energy distrbutions. one Ibr SPECT imaging and one for PET imaging, 
n is also an ot)ject of the present invention to provide a sdntillation detector pair having switchat)le triggering capabilities 
that can be used to deted and record SPECT events and also switched to perform coincidence detection for PET imag- 

so ing. It is yet anoUier objed of the present invention to employ zonal triggering capabilities for event detection in both 
SPECT or PET imaging modesw These and other ot^eds of the present invention not specifically redted will become 
dear within discussions of the present invention herdn. 

SUMMARY OF THE INVENTION 

55 

A multi-detedor head nudear camera system automatically switchade (and optimized) to perform either SPEC^ 
imaging or PET imaging is descrlt>ed. The camera system employs, in one emtxxiiment. nulti'detector configuration 
having dual head sdntillation detectors but can be implemented witti nx>re tfian two detector heads. The detectors con- 
tain switchable triggering drcuitry so that odnddence detection fbr PET imaging and non-cdnctdence detection for 
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SPECT imaging is available and can be automatically selected. Usmg a variable integration technique wifli programma- 
ble integration interval, the ev^ detection and acquisition circuitry of the camera system is switchable to detect events 
of different energy distnbution and count rate which are optimized for PET and SPECT imaging. In one mode a relatively 
larger integration interval is selected for SPECT imaging and a relatively shorter integration interval is selected for PET 

5 mode. The system also includes dual integrators on each scintillation detector channel tor collecting more than one 
event per detector at a time for either PET or SPECT mode. In PET or SPECT mode, the system also employs variable 
PMT cluster sizing having smaller di^er sizes for PET Imaging and relatively larger cluster sizes for SPECT. In PET or 
SPECT nrxxie, the system also empk^ variable centroid shape and zonal triggering. Utilizing the above programmable 
settingSi the camera system can be automatically configured to operate in either SPECT a PET imaging modes. 

10 When instructed to operate in SPECT mode, under computer control, the camera system automatically loads the 
longer integration interval into ttie detectors, programs a PMT ^Jdress tat)le to allow larger sized clusters to be gener- 
ated, queues for installation of the coilimators, and automatically s^ the trigger d^ection mode to non-coincidence. 
When instructed to operate in PET mode, under conputer control, the camera system automatically loads the shorter 
integration interval into the detectors, programs a PMT address table to only allcw smafler sized clusters to be gener- 

IS ated, queues for removal of the collimators, and automatically sets the trigger detection mode to coincidence detection. 
In SPECT mode, SPECT reconstruction procedures are applied and in PET mode, PET reconstruction procedures are 
appBed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 

Figure 1 illustrates a high level block diagram of a dual head gamma camera system employing botti SPECT and 
PET operational modes, an acquisition computer and all image processor of he present invention. 

Figure 2A1 illustrates event detection/trigger generation circuitry used by the dual head camera system of the 
present invention operating in SPECT mode of operation. 
25 Figure 2A2 illustrates event detectionAlrigger generation circuitry used by the dual head camera system operating 
in PET or coindderx^e niKxle of operation. 

Figure 2A3 illustrates circuitry of tiie coincidence timing circuit utilized by the camera system of the present inven- 
tion to provide switchable SPECT/PET triggering mechanisms for generation of valid event triggering signals. 

Figure 2B is a circuit (fiagram of the drcuitry utilized by the present invention for generating an analog total fglo- 
30 bar) energy signal. 

Figure 2C Qlustrates the preamplification/cfigitizer circuitry of the present invention (utilizing dual integrators per 
channel) for provicfing a compensated digital response signal that is digitally integrated for each PMT of tiie gamma 
detector and allowing progranmable integration intervals for PET/SPECT imaging nxxles of operation. 

Figure 2D illustrates the processing btocks of the present invention Digital Event Processor of each detector for 
35 generation of the spatial coordinates of an event local and global energy, peak PMT add^ 
output. 

Figure 3 is an illustration of a weighted centroid computation for determining the spatial coordinates of an event 
used by the present invention. 

Figure 4A and Figire 4B illustrate an overall flow diagram of aspects of the Digital Event Processor of the present 
40 invention in SPECT of PET moda 

Figure 5 is a cfiagram of a digital processor (computer system) and user interface off the present invention used, 
among other features, to automatically configure the camera system for SPECT of PET imaging nvxles. 

Figure 6A represents the light intensity response for multiple events occuning dose in time within the present 
invention. 

45 Figure 6B illustrates a f tow chart of the timing performed for dual integration on a particular channel using a pro- 
grammable integration period. 

Figure 7 illustrates an exemplary PMT array of a detector of the present invention and illustrates three (fiff^ent dus- 
ter types. 

Figure 8 is an illustration of tiie programmable PMT adcfress table drcuit of the present invention Digital Event Proc- 
50 essa. 

Figure 9 is an illustration of the x and y weight table drcuit of the present invention Digital Event Processor. 

Figure 10 illustrates a plane view of a collimator (with edges) as positioned over a PMT array of the present inven- 
tion in SPECT imaging mode. 

Figure 1 1 illustrates placement of centrally located PMTs and obscured PMTs and an exemplary strip area assod- 
55 ated witti an obscured PMT in SPECT imaging mode. 

Figure 1 2 illustrates energy absorption percentages for a typical PMT centroid anangement fbr a gamma event 

Figure 13A is a flow chart illustrating initial calibration tasks executed by the autogain procedure of tfie present 
invention as typically executed at the factory during calibration or in the field without a collimator fbr SPECT mode of 
operation. 
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Figure 13B is a flow chart illustrating routine calitxation tasks executed by the autogain procedure of the present 
invention as typicaOy executed at the site duing routine call>ration with the colGmator removed for SPECT mode of 
operation. 

Figure 13C is a flow chart illustrating routine calil)ration tasks executed by the autogain procedure of the present 
5 invention as typicaOy executed at the site during routine calbration with the collimator installed for SPECT mode of 
operation. 

Rgure 1 4A illustrates an uncompressed response of a PMT for an event verses the distance from the center of the 
event and also shows compressed responsa 

Figure 14B illustrates a summation of PMT responses (of a PMT dusteO over a one dimensional PMT configura- 

10 tion. 

Rgure 15A and Rgure 15B and Rgure 15C illustrate exenplary compression procedures realized by thedynarrac 
compression table and circuitry of the present Invention. 

Figure 1 6 illustrates an anatog voltage wave form (with k)aseline offset) of a prearrplif kation channel of the present 
Invention. 

IS Rgure 1 7A illustrates a f bw chart of the automatic t)»ellne offset calitnation of the present invention. 

Rgure 17B illustrates a fk>w chart of the software triggered baseline offset correction procedure of the present 
invention. 

Figure 17C and Figure 17D illustrate flow charts of the scintillation triggered baseline offset correction procedure 
of the present invention. 

20 Rgure 18 illustrates the memory circuit of flie present invention for provKling far PMT addresses based on an ir^ 
peak PMT address. 

Rgure 19 is an illustration of a procedure utilized by the present invention for switching between SPECT and PET 
imaging modes of operation for medk;al imaging. 

25 MULTI-HEAD NUCLEAR MEDICINE CAMERA FOR DUAL SPECT AND PET IMAGING 

DETAILED DESCRIPTION OF THE INVENTION 

In the foltowing detailed description of ttie present invention numerous specif k; detafls are set forth in order to pro- 

30 vkJe a thorough understanding of the present invention. However, it wID be obvkMis to one skilled in the art that the 
present invention may be practiced without these specific details. In other instances well known methods, procedures, 
components, and drcurts have not been desGril>ed in detail as not to unnecessarily obscure the present invention. 

Some portions of the detailed descriptions which folk)w are presented in terms of algorittims and symbolk; repre- 
sentatiors of operations on data bits wifliin a computer memory. These algorithmk: descriptions and representations 

35 are the means used by those skflled in the data processing arts to nriost effectively conv^ 

to others skilled in the art Unless specif toally stated otherwise as apparent from the fblkswing discussions, it is appre- 
ciated ttiat ttiroughout the present invention, discussk)ns utilizing terms such as "processing" or "computing" or "calcu- 
lating" or "determining" or "displaying" or tiie like, refer to the action and processes of a computer system, or similar 
electrons computing device, that manipiiates and transforms data represented as physical (electronk;) quantities 

40 within the computer system's roasters and memories into other data similarly represented as physical quantities witttin 
the computer system memories or registers or other such information storage, transmission or display devk:es. 

The various embodiments of ttie present invention descrbed herein are used In conjunction with a dual head scin- 
tillation detector camera system switchak)le t>^een SPECT and PET imaging modes of operation. With reference to 
Rgure 1, a high level diagram of a dual head detector gamma camera system of the present invention is shown. 

45 AlttK)ugh an errtxxiiment of ttie preserit invention is descrik)edwitti respect to a 

dated that the teachings of the present invention can be extended to cover systems having mae than two detectois 
(ag.. a triple head camera system is contemplated). Gently, ttie system of the present invention indudes a pair of 
gamma camera detectors 80 and 80* fdual head") conrposed of a plurality of photomultqslier tubes, PMTs, arranged in 
a two dimensional matrix convention and optically coupled to a glass plate to receive light (e.g.. visible photons) from a 

50 crystal layer 81 . The PMT array aeates a photodetector. The crystal layer can be composed of sodium iodine, Nal, and 
is typically located between a collimator (not shown) and ttie PMT array. The coffimator applied to each detector is used 
for SPECT imaging modes, kxit not for PET imaging modes, and is typically manufactured from a numt>er of holes with 
lead septas arranged in a honeyconfi> convention to collimate gamma rays ttiat strike the crystal 81 . The detector pair 
are shown in a 180 degree configuratk)n and can rotate witti respect to each other such that ttiey are at a 90 degree 

55 configuration. 

Gamma rays ttiat strike ttie Nal(ri) crystal 81 cause well known sdntillation events ttiat release a number of visible 
light photons ttiat are detected by the PMTs witti different light intensities. Each PMT reports in the form of an anak)g 
signal indicative of the amount, of light energy detected as a result of the sdntillation event. In ttie present invention, 
ttiese signals are digitized at an earty drcuit stage and are processed (figitally. The gamma camera detectors 80 and 
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80' utilized within the scope of the present invention are of the Anger type and can k>e of a number of well known and 
commercially availat)le designs arxl therefore some details of such a gamma detector will not k>e discussed in depth 
herein. However, as will be discussed herein, the present invention provides many enhancements in the image acqui- 
sition and processing of detected events so that the camera system can operate in either SPECT or PET image mode. 
5 An exemplary gamma camera detector used liy one embodiment of the present invention can contain as many as 55 
or 108 PMT& A detector of the pair can also utilize smaller diameter PMTs along the edges to increase the detector's 
field of view. An embodiment of the present invention utilizes forty-nine 76mm round PMTs and six 51 mm round PMTs 
for edge filling, however, the number of PMTs^ their sizes and their configurations can be varied within the scope of the 
present invention. 

10 It is appredated that each detector, either 80 or 80' is similarly constaicted and that discussions with respect to one 
detector are applicable to both. 

The detector pair 80 and 80' of Rgure 1 are mounted on a gantry 85 which can rotate the detectors 80 and 80' in 
various oibils around an object (patient) 1020 resting on table 87 (e.g.. tor EOT scanning operations). In either config- 
uration (180 or 90 degrees), the detector pair can rotate about a center of rotation through a number of projection 

15 angles, as is known in gamma camera technology. The gantry 85 and table 87 rest on base 89. The detector pair 80 
and 80' can also be directed transversely across the table 87 (ag.. for total body scanning operations) or placed over 
the patient 1020 for static imaging. 

Upon an event detection in either detector 80 or 80*. signals 1210 and 1212, respectively, carry trigger pulses to a 
programmable coincidence timing circuit 1050 (CTC). The CTC unit 1050 then generates valid event trigger signals 

20 over lines 1240 and 1242. respectively, for the detectors 80 and 80' depending on the mode of operation (either SPECT 
or PET). A signal carried over line 1252 indicates to the CTC unit 1 050 the proper mode of operation (SPECT or PET). 
The valid event trigger signals 1240 and 1242 are used by the detectors to start (or reset) ther accunrtulators (integra- 
tors) whKh accumulate Onte^e) the energy of the detected scintillation and are therefore called "Valid evenT trigger 
sigrials. In the PET mode, integration is not started until a coincidence is detected between detector 80 and 8a. In 

2S SPECT mode, an integration is started for each detector upon a trigger event, regardless off coincidence. After integra- 
tion and centroiding, tiie detectors 80 and 80' output an X. Y. and Z value over lines 1220 and 1222. respectively. These 
signals indicate the coordinate vahie (X. Y) (localization') of the detected event with respect to the detector .and its 
measured energy value, Z. 

WHhIn embodiments that utilize more than two detector heads, event detection infonnation from each detector 
30 head is fbnwaided to the CTC unit 1050 that then detects ooincUence between any twodetectors feeding the CTC event 
detection information (when in PET imaging mode). In SPECT mode, each detector reports event information in non- 
coincidence in a similar ^ion as the dual detector system. 

AHhough the location of such hardware is not material to the present invention, hardware is included wittiin each 
detector for di^zing the PMT channel signals, correcting them^ and outputting X. Y, and Z values as will be discussed 
3s further betow. To this extent each detector 80 and 80' contains preairpfific^ 

Digital Event Processor. These are discussed further below. This hardware can be tocated within or outskie the sdntit- 
lation detector 80 and 80*. 

The values transmitted over bus 1220 and 1222 are input to an acquisition computer system 1055 whu:h contains 
a general purpose digital computer system 1112 (Rgure 5) as descrbed to fbltow. The acquisition computer 1055 

40 stores the values for each detected event for each pn^ection angle and this information is routed to an image processor 
1060 whwh contains a standard user interface. The user interiaoe provMes a user input devtee for indk»ting which 
mode of operation (e.g. SPECT or PET) is requested. This device can also be kxsHed in computer 1055. The image 
processor 1060 performs image reconstruction and attenuation and uniformity correction based on the acquired event 
data for SPECT or PET imaging modes. The image processor 1 060 is also coupled to a display unit 1065 (whk:h can 

45 include a hardcopy device) for visualizing images captured by the camera systera 

DFTgCTOR ELECTRONICS 

The circuitry and logic of each detector 80 and 80' for signal processing is further illustrated within the following f ig- 
50 ures: Rgure 2A1 . Figure 2A2. Rgure 2B, Rgure 2C and Rgure 2D. It is appreciated that this circuitry can alternatively 
be k)cated outskie the detector heads 80 and 80' or partially spread inskJe and partially outskte of the detector heads 
80 and 80'. It is further appreciated tiiat the acquisition computer 1055 also comprises a general purpose computer sys- 
tem 1 1 12 (see Rgure 5) as will be discussed in more detail. An embodiment of the present invention can also be imple- 
mented utilizing disaete electronic components in lieu of a general purpose conputer system. 

55 

TRIGGER DETECTION AND SWITCHABLE PET AND SPECT EVEffyrT RECOGNITION 

Since the present invention system is switchable between SPECT and PET imaging, ttie system employs two sef>- 
arate trigger mechanisms for recognizing a valkl event and generating a valkJ event trigger signal. There are at least 
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two dfferences between these trigger mechanisms for the SPECT and PET imaging systems. Generally, the first differ- 
ence is that the actual trigger detection circuitry is different for SPECT and PET modes of operation t>ecause the count 
rate is faster for the PET system. The SPECT mode utilizes a leacfing edge discriminator in its event detection drcuitry 
while the PET mode utilize a constant fraction discriminator in its event detection circuitry t>ecause the time t>etween 
5 valid events from the two detectors needs to be known with high accuracy. Secondly, the PET mode utiGzes a coinci- 
dence circuit to determine a valid event from a detected event while tfie SPECT mode does not These drfferences will 
be explained below. 

Figure 2A1 illustrates detection logic 100 for detector 80 utilized in the SPECT mode. A duplicate circuit of 100 
(lower) is used for detector 80' fbr SPECT mode. In a system having more than two detector heads, logic 100 is dupli- 

10 cated fbr each detector head. The outputs 1 130 and 1 132 from these two circuits belong to buses 1210 and 1212, 
respectively, and are forwarded to the CTC unit 1050 (see Rgure 1). Regarding Figure 2A1 , the trigger signal genera- 
tion circuit 100 is diplicated for each detector as shewn. Therefore, discussions regarding the operation of circuit 100 
wHh respect to detector 80 are applicable to the duplicate circuit for detector 80' of the present invention. According to 
the arrangement of the PMTs in tfie detector, the output of each PMT associated with a particular spatial quadrant (or 

16 "zone") of the detector rnatrix is sent to one of four trigger detection drcuHs 110a 110b, llOcor llOd. Eachzonerep- 
resenls a separate spatial portion of the scintillation detector. TTie zones can be divided horizontally aaoss the imaging 
suriace. vertically, or both. Although the precise alignment of each PMT in each zone is not critical to the present inven- 
tion, the zones can be overiapping. The PMT signals shown A0-A54 are voltage signals. The signals of the first zone, 
AO... A17 are sent to trigger circuit ("event detection") 1 10a where each is coupled to a delay line 102 and an amplifier 

20 1 04 which together create a 200ns dip drcuit The trigger pulse received from the PMT drcuits is roughly a 200 ns ana- 
log pulse. An event occurring within a particular zone will cause an event detection drcuit of the corresponding zone to 
generate an output signal or "event indcation." In ^steins having more than two detectorSi the above event detection 
drcuitry is located in each detector. 

For the SPECT mode of operation, the dipped signal from circuit 1 04 is ooi4)led to the positive end of a discrirn- 

25 nator drcuit 106 and a computer controlled r^erence input is coupled to a threshold input droflt 108. The reference 
signal at 108 is coupled to receive the output of a corrputer controlled DAC (not shown). Therefbre. only trigger signals 
over the threshold voltage are allowed to pass through drcuit 1 06 and they are dipped to 200ns. The output of tiie com- 
parator 106 is tfien coupled to the input of OR gate 120 via line 122. Line 130 will assert a triggering pulse when ever 
a PMT of the designated zone detects an event This drcuitry is replicated for each of the other four zones of the detec- 

30 tor PMT matrix for detector 80 (e.g., signals All . . . A30 feed circuit 1 1 0b which generates a trigger over line 1 24, sig- 
nals A24. . .A43 feed drcuit 110c which generates a trigger over line 126. and signals A37 . . . A54 feed drcuit llOd 
which generates a trigger over 128). The same is true for detector 80*. It is appredated that additional or fewer event 
detection drcuits (e.g., for more or fewer zones) can t>e used within the scope of the present invention. Multiple trigger 
channels (zonal triggering) are used to maintain a Ngh signal to noise ratio in the presence of correlated noise arxl to 

35 maintain low dead time at high count ratea The trigger channels are cveriapped, as discussed above, to prcR/ent sen- 
sitivity loss at the zone boundaries. 

Trigger lines 122, 124. 126. and 128 of Figure 2A1 are coupled to the input of OR gate 120. Therefore, when an 
event is detected by the camera detector 80. OR gate 120 of trigger drcuit 1 10 of the present invention will generate a 
trigger pulse over line 1130. These triggering pulses are used, along with otfier signals, to generate Start(t0) and 

40 Start(t1) and (after going through CTC drcuit 1050) are used the integration drcuits 280(0)-280(54) (Rgwe 2C) of 
the present invention to start integration of the PMT signals for the detected ganvna event. When an event is detected 
t>y the camera detector 80*. duplicate OR gate 1 20 of the duplicate trigger drcuit fbr detector 80* of the present invention 
will generate a trigger pulse over line 1 132. Lines 1 132 (fbr d^ector 80^ and 1 130 (fbr detector 80) eire within buses 
1212 and 1210. respectively (see Rgure 1). 

45 Event detection fbr the PET imaging mode comprises different electronics which are illustrated in Figure 2A2. Cir- 
cuit 100' is associated with detector 80 and aduplicate drcuit (lower) is associated with detector 80*. In a system having 
more than two detector heads, logic 1 00' is duplicated for each detector head. Unless discussed as different, drcuit 1 00' 
(used for PET) is analogous to drcuit 100 (used for SPECT). The outputs 1130* and 1 132* from these PET drcuits 100' 
belong to buses 1210 and 1212. respectively, and are fonivarded to the CTC unit 1050 (see Figure 1). The ev^ detec- 

50 tion logic 1 1 0a - 1 0Od* for each detector zone utilizes a constant fraction discriminator when in PET mode in lieu of the 
leading edge discrirrvnator fbr SPECT imaging used in drcuit 100. This is done to rraintain timing accuracy when per- 
forming PET imaging. An event occurring within a particular zone will cause an esrent detection drcuit (1 1 0a - 1 lOd*) of 
the corresponcfing zone to generate an output signal or "event indication. " The output of the event detection drcuits 
1 1 0a'-1 OOd' are dpped to one-half the width of the desired coincidence timing window. For each detector, the outputs 

55 Of the event detection drcuitry are fon^varded to an OR gates 1 20b. The outputs 1 1 30* and 1 1 32* of these gates are for- 
warded over buses 1210and 1212, respectively. It is appreciated that pipelining of the preampli^ 
(as will be discussed below) allows a 240 ns "memory" therefore aDowing capture off the entire pulse that generated the 
detected coincidence pulse in PET imaging mode. 

Figure 2A3 illustrates the logic of the CTC unit 1050 (as shown in Figure 1). Bus 1210 carries signals 1130 and 
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1 130' while bus 1212 cairies signals 1 132 and 1 132*. Signals 1130' and 1 132' are input to AND gate 1 125 which gen- 
erates a signal ever line 11 35 only if signals 1 130' and 1 132' are in coincidence (ag., within a 1 0-1 5 ns window). Signal 
1 1 35 is coupled to multiplexer 1310 and to multiplexer 1320. Sigial line 1 130 is coupled to multiplexer 1310 while signal 
1 1 32 is coupled to multiplexer 1320. A control (or select) signal 1252 is coupled to the select inputs of multiplexers 1310 

5 and 1320. When the control signal 1252 selects PET imaging operation mode, then the signal over line 1 135 passes 
over both line 1 240 to detector 80 and to line 1 242 to detector 80'. These are the vafid event trigger signals. When the 
control signal 1252 selects SPECT mode, then the signal over line 1 130 is carried over Kne 1 240 to detector 80 and the 
signal over line 1 132 is canied ever line 1242 to detector 80'. These are the valid event trigger signals. It is appreciated 
that the signals over line 1240 are used to trigger the integrators associated with the preamplrfication and digitization 

10 circuitry for detector 80 while the signals over line 1 242 are used to trigger the integrators associated with the pream- 
plification and digitization circuitry for detector 80'. The signals carried over line 1240 aixJ 1242 represent "Valid* event 
detectioa The nrultiplexers 1310 and 1320 allow a switchaksle mode of operation t>etwe^ event detection for SPECT 
and PET imaging modes under the present invention including a coincidence detection mode (line 1135) and a norv 
ooincidenoe detection mode (lines 1 130 and 1 132). 

IS For detector 80, the valid event trigger signal 1240 is used as trigger signal 130 to convey start (tO) and slart (t1) 
while for detector 80' the valid event trigger signal 1242 is used as trigger signal 130. 

In a system having more than two detector he^, in SPECT mode, Itte event indication signals from each detector 
are routed back to their associated detector as the valid event trigger. In PET nrKXie, a series of AND gates are coupled 
so that each possible connbination of two detecta is represented by an AND gate and, in such way, a coincidence 

20 between any two detectors of the system can be determined. Once a coincidence is detected, a valid event trigger is 
forwarded t>ack to the two detectons in ooinddenoa 

ANALOG SUM OF GLOBAL ENERGY 

25 Figure 2B illustrates an analog summing circuit utiized l)y one embodiment of the present invention. This circuit is 
included in detector 80 and a separate drcut in detector 80'. The voltage signals AO. . . A54 from each PMT channel 
(e.g., signal from each PMT) are summed together and output t>y amplifier 142. An offset voltage is fed into the sum- 
ming circuit via input 144. The offset voHage is controlled by a computer controlled DAC. The output, or total energy of 
the gamma event, is generated over line A55. This output represents the analog sum of all of the voltage si^ials from 

30 each PMT channel plus an adjustable offset from circuit 144. The analog signal over ASS is called the analog global 
energy signal. This glok)al energy signal may k>e supplied to an available channel (ag., channel 55) for preamplif ication, 
digitization and integration by drcurt 200 (Rgure 2C). The Digital Event Processor of each detector the present inven- 
tion receives signal A55 of its corresponding detector. As will t>e discussed to follow, the glot>al energy value is also 
computed digitally by summing the digitized integrated channel signals of the PMTs. Either of these methods can be 

ss utilized within the scope of the present invention. 

PPEAMPUFICATIQN PtemZER 

Refer to Figure 2C which illustrates the prearrplif ication digitizer circuits 200 of the present invention for each of the 

40 55channels(plusonefor ASS) Ibr detector 80. A dijplicate of this drcuit is used wi^ 

system having more than two detectors, circuit 200 is duplicated for each detector. Circuit 200 is used for both SPECT 
and PET imaging mode and utilizes switchable circuitry to optimize t>etween modes. These circuits perform the pream- 
plif ication, digitization, and integration for each analog voltage signal for a PMT channel . The output of circuit 200 is fed 
to a Digital Event Processor for detector 80. A duplicate preamplffication digitizer circuit 200 is used for detector 80* and 

45 outputs to a duplicate Digital Event Processor of detector 80*. As will k>e discussed, each preamplif ication digitizer circuit 
for each channel contains two separate integrator circuits. Unless discussed othenwise, circuit 200 for detector 80 is 
analogous in operation and structure to its duplicate circuit of detector 80*. Circuit 280(0) corresponds to the current 
(amps) output signal received directly from PMT #0 (ag., channel 0) and this circuit 280(0) is separately replicated for 
each of the 55 PMT channels of the present invention and, as shown, circuits 280(0) to 280(54) operate to simuHane- 

50 ousiy process the current output signals for PMTO to PMT54. Regarding circuit 280(0), tiie current signal output of 
PMTO is fed into a current to voltage converter 210 and the output of this signal is fed through a resistor to a voltage 
gain amplifier 222. 

A computer controlled digital to analog converter (DAC) outputs two adjustment signals over line 212 and 214 for 
k)aseline voHage oorrectioa The signal over line 212 is fed through a resistor 21 6 for coarse adjustment and line 21 4 is 
55 fed through resistor 21 8, which has much larger resistance (e.g., on the order of 200X) than resistor 21 6, for fine adjust- 
ment. The signals received via lines 212 and 214 provide a baseline offset voltage adSustment to the output signal 
received from ttie PMTO. A conrputer controlled digital to analog converter (DAC) outputs a voltage adjustment signal 
to drcuit 220 to control the gain of amplifier 222 fiaving an exenplary gain adjustment of 10:1 . The analog gain adjust- 
ments are coarse acQustment with fine gain adjustnrients performed by the calibration table, see Figure 2D as discussed 
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to fdlm. Ihe baseline offset adjusted signal and Hie gain adjusted signal Is tfien output at point 224 as signal AO for 
each channel. Similarly, for each PMT channel, the atxve circuitry is replicated for generating signals A1 to A54. The 
trigger signal 130 Is uniformly supplied to each of the circuits 280(0) to 280(54) so that each channel is triggered coin- 
ddentally. 

The output of gain amplifier 222 is then fed into resistor 226. The voltage input 228 and voltage input 230 are cou- 
pled tfrough respective resistors to the output of resistor 226. The output of resistor 226 is then fed into amplifier 232 
and capacitor 234 in series. The outputs of the amplifier 232 and capacitor 234 are then coupled to the input of the ana- 
log to digital converter (ADC) 236. The above circuit (ag., from point 224 to the input of ADC 236) is utilized for pulse 
insertion used for diagnostic purposes that are not particulariy pertment to the present invention. Pulses may be artifi- 
cially inserted via inputs 228 and 230. 

Referring to Figure 2C. the ADC 236 converts the analog signal AO to cfigrtal samples based on the frequency of 
dock input as shown. One embodiment of the present invention utilizes a sample frequency of 20 to 25 MHz as a sam- 
ple dock. The output of ADC 236 is then fed to the input d tm adders fimegrators" or "a^^ 
coupled in an integration configuration. The present invention utilizes dual digital integrators for each PMT channel in 
order to more effectively process concfitions wherein two gamma events are detected in dose temporal proximity. Each 
integrator 238 and 240 (of each channel 0-54) contains a separate register (accunulator) for containing the current sum 
value and perfomns piecewlse linear integration with a dip binary value at 1023 (ag., no roDover is allowed). The out- 
puts of botii integrators are coupled to multiplexer 241 . The multiplexer 241 selects between one of the two integrator 
registers (accumulators) for output to the latch drcuHs 242 and 244. Latch drcuits 242 and 244 comprise a two stage 
FIFO anangement. The trigger pulses Start(tO) and Start(t1) received over 130 are used to reset the integrators. These 
trigger pulses arrive over line 1240 for detector 80 and over line 1242 for detector 80*. At the end of a programmable 
integration (accumulation) period, the value of the integration process for eHher inte^ator is ttien stored wittiin a two 
stage latch drcuit of 242 or 244. The output of the latch droit 242 or 244 is the digitized value of the signal generated 
by P^m1br a spvenevem and ttiisvafoe is designated as DO. The (j^^ Digital Event 

Processor 300. The above is explained in more detail with reference to Figure 6B. 

Under computer control, the duration of the integration period for all integrators is programmable over line 235 
which is coifHed to a computer processor. The programmed period (having a resolution of the sampling dock frequency 
of 20-25 MHz) is loaded into register 1 330 over line 235. Upon a valid event detection signal (over fine 1 240) that resets 
integrator 238, a corresponding counter 1332 is also reset Counter 1 332 is coupled to the sampling rate dock. The out- 
put of the counter 1 332 is coupled to a comparator 1 335 which also receives an output from register 1 330 which con- 
tains the programmable integration (accumulation) period. Upon a coincidence, an output signal is can-led over line 
1 340 which signals the completion of the integration period for integrator 238. Upon a valid event detection signal (over 
line 1240) that resets integrator 240, a con-esponding counter 1337 is also reset. Counter 1337 is coupled to the sam- 
pling rate dock. The output of the counter 1 337 is coupled to a comparator 1 342 which also receives an output from 
register 1 330 wNch contains the programmable Integration (accumulation) period. Upon a coincidence, an output sig- 
nal is carried over line 1341 which signals the completion of the integration period for Integrator 240. H is appreciated 
that like the integrator drcuits 238 and 240, the counters 1332 and 1337 are not effected tiy the trigger pulse unless 
tiiey are kfle at the time the valid event trigger pulse is received. 

It is appreciated that register 1330, counters 1332 and 1337, and comparators 1335 and 1342 are coupled to each 
channel drcuitry 280(0)-280(54) and the connections shown with respect to channel drcuitry 280(0) are exemplary of 
the connections made for each other channel. In other words, only one register 1 330 is utilized to contain the integration 
interval associated with all channel integrators for a given detector. Specfffoally, the output of comparators 1335 and 
1342 are coupled to the accumulators 238 and 340 of all channel drcuitry 280(0)-280(54) in the fashion shown in Figure 
2C with respect to channel drcuit 280(0). 

The integration (accumulation) period programmed into register 1330 is programmable with a resolution of 20-25 
MHz in one embodiment of tiie present invention. Durhg PET imaging, the accumulation period is set to approdmately 
320 ns per event due to the high count rate received during PET imaging and the high energy distribution of ganrvna 
events that are detected (51 1 K electron-volts) in PET imaging. During SPECT operational modes, the accumiJation 
period is set to approximately 840 ns as a result of the Icwer count rate detected and the fower energy distribution of 
gamma events detected. Therefore, using the above drcuitry, the present invention provides a programmak)le mecha- 
nism for optimizing the preamplif ication/digitizer drciflts for eittier PET or SPECT modes and this mechanism is switch- 
able. It is appredated that a dupticate drcuit 200 is induded for detector 80' but its valid event trigger signals (ag., 130) 
arrive over line 1242 instead of line 1240. The program signal 235 is coupled to both drcuits 200 such that both detec- 
tors 80 and 80' are simultaneously configured for PET or SPECT mode of operation. 

Using the dual integrators of Figure 2C of the present invention, either one or both of the accunulators of 238 or 
240 can be enabled to integrate an incoming signal from the output of ADC 236. The two integration results are multi- 
plexed onto a common data patti and eittier result can be selected and stored in the two stage latch drcuit. Each inte- 
grator can be separately triggered by valid event trigger 8tart(tO) and 8tart(t1). In operation, when a trigger signal 
occurs, if either integrator is available (ag.. not integrating and not hokfing an intefi^ed result) then that accumulator 
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is reset and enabled to t)egin integrating the e/ent. When either of the accumulators comptetes, the integrated value is 
transfen-ed to the first FIFO stage (e.g.. latch 242) assuming this stage is availat)le. If it is not available, the accumulator 
holds the value. Integration continues for a predetermined period of time after the trigger signal until a sufficient amount 
of the gamma event's energy is integrated. Values are transferred from FIFO stage 1 (latch 242) to FIFO stage 2 (ag., 

5 latch 244) as FIFO stage 2 becomes available (e.g., transfer its value). When data is written to FIFO stage 2, the 
present invention signals that data is ready to be transfenred to an associated Digital Event Processor (DEP) 300. see 
Figure 2D. Each detector 80 and 80* contains its own DEP circuit 300. 

Refem'ng to Figure 2C, the dual accumulator design of the present invention provides an implicit mechanism for 
hancRing event pile-t^) where two events interact during the same time period. Since the present invention computes 

10 positions based on local PMT dusters. pile-i|> events which occur in different regions of the detector can be property 
positioned and such are called temporal pile-ups. If the two events that are involved in a temporal pilenip happen to be 
separated k)y more than the trigger channel deadtime. then both accumulators will be enabled and both events vm'll k>e 
fully ffitegrated. Accuracy of positions will be impacted by the spatial distance separating the temporal events. The 
neater the separatioa the lower the impact. This is discussed in more depth below. 

IS Within circuit 200. the circuit 280(0) is replicated for each channel as shown in Rgure 2C. The output currerrt sig- 
nals from PMT#1 through PMT#54 are fed into circuits 280(1) through 280(54). The digital data signals DO to D54 are 
output from circuits 280(0) through 280(54). respectively. Each of the 55 preampllfication digitizer drcuits 280(0)- 
280(54) are coupled to receive valid event trigger signals start(tO) and start(t1) from line 130. It is appreciated that an 
extra channel (e.g.. a preamplification digitizer drcuit 280(55)) can be added in order to process the analog global 

20 energy signal from ASS (see Rgure 2B). In this embodiment, the output D55 would correspond to tiie amplified, digi- 
tized and integrated value for all channels (ag., ttie digitized value of the analog glok>al energy of tiie event). In such an 
embodiment, the value D55 would be output to the DEP 300 (as vifill be discussed to follow) with an appropriate PMT 
address value indicating the data as analog global energy data. 

The preamplification drcuits of Rgure 20 can t>e directly adjusted using gain (ag., line 220) and baseline offset 

25 (e.g., lines 212 and 214) adjustments. The atxive atjQustment lines are referred to as control link signals which are cou- 
pled to computer controlled addressable DACs. The amourrt of nominal baseline and the amount of variation associated 
with levels of adjustment can be i^ed to determine the accuracy of. given channel. 

At the completion of the programmed integration period, when latch 244 data is present, then all of the digital data 
stored in each second stage latch for each channel is transferred to the appropriate DEP 300 (Rgure 2D) over bus 307. 

30 This is caUed a data "Iransfer" to the DEP 300 and occura for bolh detectors 80 and 80*. 

DIGITAL EVENT PROCESSOR 

Refer to Rgure 2D which illustrates drcuitry 300 of the DigrtEril Event Processor (DEP) of the present invention. It is 

35 appreciated that the digital processor 11 12 is ricl part of the DEP 300 drcuit t)ut is rather 

computer system 1055. Discussions herein de8crit>e tiie DEP 300 assodated with detector 80. A duplicate DEP drcuit 
300 assodated witti detector 80* operates in an analogous fashion with respect to output it receives from ttie preampli- 
fication/digitizer drcuits associated with detector 80*. The duplicate DEP 300 assodated with detector 80* is also simi- 
larly coupled to the digital processor 1112. The output of DEP 300 for detector 80 is carried over bus 1220 to the 

40 acquisition corrputer 1055. The output of duplicate DEP 300 for detector 

computer 1 055. Each DEP 300 of the detector pair 80 and 80' operates IrKlependenfly of the ottier to generate the X, 
Y, and Z outputs to the acquisition computer system 1055. 

The following discussion descrit)es the DEP 300 assodated with detector 80. however, it is applicable to the dupli- 
cate DEP witttin detector 80*. The digitized and integrated signal values for the PMT channel over Ones DO through 

45 D54 are fed over bus 307 to FIFO 310 and to calibration table 315 shown in Rgure 2D. The data over bus 307 repre- 
sents the digitized integrated signals supplied from each PMT channel drcuit, 280(1), in response to a triggering event 
such as a sdntillation event. The dgitized data over bus 307 is stored into a raw view FIFO 310 that can be accessed 
over bus 397 by the digital processor computer 1112. The raw view FIFO allows data to be pulled from tiie input data 
stream wittiout intenupting the normal data flow from the PMTa This is utilized for on-the-f ly baseline acjOustment A cal- 

50 ibration table 315 receives as inputs (1) the digitized integrated channel data from bus 307 and also receives (2) PMT 
address (ag.. indicator) mtvtoer over bus 302 of the reporting PMTs in order to correlate the digital data over txjs 307 
with tiie proper PMT channel output from drcuits 280(i). 

The calibration table 315 contains a lookup table (LUT) for providing a gain output which can vary witti the PMT 
number input (represented in one embodiment as an address); this spatially variant gain is applied to the integrated sig- 

55 nal value received over bus 307 and tiie result is output over bus 347 which is a corrected or finely corrpensated inte- 
grated signal value for each PMT channel. The gain value stored in the calibration table 315 is a fine gain adjustment 
dependent on the PMT number whereas ttie gain amplifier 222 is a coarse gain aci^ustment The calibration table 315 
also provides a baseline adjustment computation by digitally subtracting ttie reference voltage inserted tiy the baseline 
offset drcuitry (e.g., inputs 212 and 214) of ttie preamplification drcuts for each PMT channel independenfly. The out- 
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put over bus 347 is the (figitized values supplied from 307 with this baseline adjustment. 

The output of the calbration table 315 is transfenred over bus 347 to a peak detect circuit 320 which analyzes all 
the calibrated results of all 55 channels (for a given data transfer) and selects the PMT number having the largest inte- 
grated channel signal (ag., "energy^ output fbr a given measured event; this is the "jpeak PMT* The maximum inte- 

5 grated signal value and the associated channel address are retained for later use m the DEP process of the present 
invention. The Integrated signal associated with the peak PMT is output from peak detect circuit 320 over bus 31 7 as 
value PD. The PMT address number associated with the peak PMT is output by circuit 320 over bus 31 2 as value PA. 
To support anafog global energy data being transfen'ed over a digitizer channel of 280(i), the peak detect 320 can be 
disabled for a given PMT address. Bus 347 is also ooi4)led to a gfobal energy accunrulation (GE accumulation) circurt 

10 330. Circuit 330 sums the corrected integrated channel signals output over bus 347 for each of the PMT cfiannels for a 
given event. The output of circurt 330 is the digitized global energy GE (which is a digital sum of all of the PMTs digital 
irtegrated signals) and is transterred over bus 322 whk:h fomfis output GE and ato 
sion table 355. 

Buffer 325 of Figure 2D stores the digital Integrated signal vafoe of each PMT eta 

IS ate PMT address received over bus 347. Buffer 325 can be irrplemented in RAM or other memory storage devk:a The 
PMT integrated signal values for all channels are stored in buffer circuit 325. The address of the peak PMT is output 
over bus 31 2 to drcuit 335 which is the PMT address taUa Circuit 335 contains a lookup table that outputs a PMT dus- 
ter based on a peak PMT address input from bus 31 2. The PMT duster is a collection of PMTs whose integrated chan- 
nel responses (as well as a total energy fbr the event) are used to perform the DEP computation to determh^ 

20 location of the event (e.g.. the centrokJ of the duster). By using a fookup table at drcuit 335. the present invention is able 
to provide spatially variant duster shape and to vary spatially the number of PMTs tiiat form a given PMT duster. When 
operating in PET imaging mode, drcuit 335 is programmed to output relatively smaller (ag. 7 PMTs) duster sizes wNle 
operating in SPECT imaging mode, drcuit 335 is programmed to output relatively larger (ag. 17-19) sized dusters. 
WHhin the scope of the presem invention, fbr an input peak PMT address, the shape off the res^ 

25 and the number of PMTs that make up the selected PMT duster vary based on the spatial tocation of thepeak PMT 
address witNn the cverall PMT matrix. It is appreciated that in SPECT imaging mode the PMT address table 335 can 
vary the PMT duster associated with the peak PMT address based on a selection for high or ksw resolution wherein a 
low resolution mode may require 7 PMTs per PMT duster where a higher resolution nfKXie may require 9 to 19 PMTs 
per PMT duster. Therefore, a resolution indication signal (not shewn) is also input to the table 335. It is appredated that 

30 in lieu of a resolution incficator signal, the entire table 335 can be reloaded witii data fbr cfifferent desired resolutions. In 
such case the resolution signal is not used as an addressing signal but only initiates the downloading off the new infor- 
mation. 

A sequence counter 390 is coupled to ttie PMT address table 335 via bus 367. The PMT address table 335 controls 
(1) the nurrfoer off PMTs in the selected PMT duster fbr the given evem and (2) the type valu 

3s is then stored in circuit 340 and hekJ tturoughout the spatml computation for a given event). The PMT address table 335 
also contains the address of the analog gfobal energy channel. The sequence counter 390 then counts, sequentially, 
from one to the numt>er of PMTs that are associated with the selected PMT duster and sequentially presents each 
count value over bus 367. In one embodiment, the PMT address tal3le 335 is HseH addressed by two values, (1) the 
MSB of tfie address that originates from the peak PMT address value over bus 312 and (2) the L^B of the address value 

40 thai originates from the ooum value of the sequence counter 390 over bus 367. The last entry ^ 
table fbr a given peak PMT includes a stop code thai indicates the end of the PMT duster 

address. The centroid computation drcuitry ttierefbre stops (e.g., is terminated) when ttie stop code is reached (or 
equlvalentiy when the maximum count value is reached as reported tyy the table 335). 

The PMT address table 335 outputs ever bus 372, in sequence based on the sequence counter 390, the PMT 

45 addresses of each PMT of the PMT duster used in the spatial oomputatfon for a given everiL The order in wh^ 

PMTs are presented over bus 372 is gcvemed tyy the fookup table stored in the PMT address table 335 based on the 
peak PMT address value from bus 312 arxl the count value over bus 367. The PMT address values output from drcuit 
335 are also coupled to data buffer 325 to address this memory drcuit which will output the appropriate integrated chan- 
nel signal value over bus 352 for tfiat PMT. This output is used in the DEP*s spatial computation. 

so Referring still to Figu'e 2D, the PMT cluster type value is output from the PMT address table 335 to the memory 
drcuit 340 which hoUs the PMT duster type value throughout a centroid (e.g., coordinate) computation. The PMT 
address values over bus 372 are also fed to a weight table drcuit 345 that contains a lookup table correlating PMT 
address value with a given x and y weight value used for the spatial corrputation drcuitry of the present invention. By 
providing the lookup table within drcuit 345. the weight assodated with a given PMT depends on its address value and 

55 is oonrelated with the spatial location off the PMT and it is also dependent on the PMT duster type value from bus 377. 
The weight value assodated with a given PMT address can also vary t>ased on the type vafoe stored in ^ 
is based on the peak PMT for the PMT duster of the given event; this is accomplished by the type regstration drcuit 
340 off the present invention. Therefore, if tiie peak PMT address corresponds to be an edge or corner PMT. then the 
PMT duster wOl be off a spedal type and the weights assodated wifli the PMT addr 
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can be adjusted to aocouit for the missing PMTs off the centroid computation (e.g., the PMTs that are not available due 
to the edge or corner location of the peak PMT). Therefore, the PMT waght value output from circuit 345 depends on 
(1) the PMT address value and (2) the type value from circuit 340 that is based on the address value of the peak PMT 
for the PMT cluster. The PMT address table 335 defines the PMT addresses that constitute the selected PMT cluster. 

5 In one embodiment the type registration circurt 340 may contain a k>okup table based on the peak PMT address 
value output from drcuit 335. As discussed, the PMT cluster type value is an offset into the weig^ table 345 for a given 
PMT cluster and is used to pfiovkie variations in the weights output of the weight table 345 for a given PMT address 
value of a given PMT cluster. Circuit 340 is coupled to supply the offset value to drcuit 345 via bus 377. According to 
the operatk)n of the present invention, for a given PMT cluster, the PMT addresses that make ip the PMT cluster are 

10 sequentially output over bus 372 to the weight table and a constant type value is generated and output over bus 377 
(based on the address of the peak PMT). The weight table 345 then outputs an x weight over bus 362 and a y weight 
over bus 357. for each PMT of the given PMT duster. 

The X weight value over bus 362 of Rgure 2D supplies a multiplier accumulator drcuit (MACx) 370 for the x coor- 
(finate computation and the y weight value over bus 357 supplies a multiplier accumulator drcuit (MACy) 365. These 

IS drcuits 370 and 365 are reset and initialized at the start of the spatial conriputationfr^ 

weight values are utifized in the spatial computations for a given event (ag., gamma interaction) and specify tfie amount 
of contributfon a given PMTs integrated channel signal value should canry (for a given PMT of the PMT duster) in the 
coordinate computatfon process. 

As discussed above^ the PMT addresses of the PMTs off a given PMT cluster are placed over bus sequence. 

20 Bus 372 of Rgure 2D is coupled to address the buffer 325 with the PMT address value for each PMT involved in the 
PMT cluster. In response to the address value of a given PMT of the PMT duster, buffer 325 outputs the stored and 
corrected digital signal value assodated with the PMT (for a given event) as received over bus 347. This conreded sig- 
nal value is transferred over bus 352 to btock 350 to block 355 and to subtractor 391 . A 
the digital signal vahjes of the PMT^ involved in the given PMT cluster a^ 

2S energy (LE) value which is generated over bus 337. 

Btock 350 is a buffer for containing either a digital or analog global energy value when operating in a mode wherein 
the global energy data is transferred over a digitizer channel and stored in the storage buffer 325. In this mode the glo- 
bal energy stored in btock 350 is a digitized version of the analog gfobal energy signal. Output from the analog signal 
ASS (of Figure 2B) Is fed to a preampliftoatfon (figitizer channel off 280(i) and then over bus 307 to the calbration table 

30 315 and staed in buffer 325. The gfobal energy nray also be computed by accurnulator 330 by 

integrated signals from each PMT channel. The gfobal energy value is then stored in buffer 350. Therefore, the gfobal 
energy value GE over fine 322 can be (1) a dgitized value of the analog gfobal energy value or (2) a digital summed 
value of tfie digital signals of each PMT channel. 

The dynamfo oompressfon table 355 of Rgure 2D receives the global energy value of the detected gamma event 

35 over bus 322 and also receives the digitized integrated channel signal value for a given PMT of the PMT duster ever 
bus 352. The dynamic compression table 355 contains a fookup table off compensalfon values for the digitized inte- 
grated channel signals. The output of the fookup table is driven on bus 392 to a sut)tractor 391 which also receives the 
signal value over bus 352. The table 355 receives as an address the MSBs of the global energy (from bus 322) and 
predetermined bits of the signal value for each PMT over bus 352. Via the connection with the subfractor 391 . the signal 

40 data over fine 352 is sNfted left by four bits (ag., multiplied by 16). The output from table 355, in one embodiment is 
subtracted from this left shifted signal valua The output of the subtraclor drcuit 391 is the dynamfo conpressed inte- 
grated signal value for a given PMT channel and is then fed to the MACx 370 drcuit, the MACy drcuit 365 and also to 
an energy multiplier accumulator (MACz) drcuit 360. 

The dynamic conpression tat)le 355 of Rgure 2D is utilized by the present inventfon to alter the integrated channel 

45 signal output from a PMT channel so that when summed with altered signals from other channels, the summatfon signal 
will be more linear in nature. The infbrmatfon stored in the dynamfo compression table 355 that is used to perfrym the 
signal conversion is readily programmable witNn the present invention and different conversfon data sets may be stored 
(down foaded) in the table 355 at one time. Since the conversion data may be altered readily (e.g.. a new set can readily 
be downloaded, if needed) the dynamfo oompressfon table 355 of the present invention is nrxxlif iable. The subtraction 

50 logic391 is a part of the coripresdon procedure used tyy the present inventfon and is used in or^ 

ory size requirement of the compression table 355. Therefore, it is appreciated that given a larger memory size, the 
subtracter 391 may be eliminated from the present inventfon and integrated into the memory 355 by altering the data 
stored therein. The output of the dynamic compression table 355 is called the 'dynamic corrpressed" or "compressed" 
integrated signal data lor a partfoular channel and is supplied to the cenfroid computation fogic via bus 387. 

55 Therefore, as the sequence counter 390 counts through the PMTs of the PMT duster, the buffer 325 supplies the 
integrated signal assodated with each PMT. Also, the weight table 345 supplies the x and y weight values associated 
with each sequenced PMT The MACx drcuit 370 multiplies the x weight value and the dynamfo compressed signal for 
each PMT and accuirulates these values for each PMT of the PMT duster as the sequencer counts. The MACy drcuit 
365 nudtipfies the y weight value and the dynamic cornpressed dgnal for eac^ 
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each PMT of the PMT cluster as the sequencer counts. The MACz circuit h^ two mputs. one is coi4)led to bus 342 
which is programmed to a value of "1 r and the other Input Is coupled to t)us 387 and therefore will accunrulate the inte- 
grated signal of each of the PMTs of the PMT cluster to generate a value of the local energy (LE) over bus 337. 
After the sequencer 390 of Figure 2D reaches the last PMT of the PMT cluster, the MACz circuit will output the com- 

5 plete LE value over bus 337 which is coupled to a 1/LE drcuit 385. This is a lookip table that provides the inverse (e.g.. 
(LE)'^) of the LE value. Drcuit 385 can also be reafized using a divider drcuit. The value of (1/LE) is then output over 
bus 382 to X multiplier circuit CMULx) 375 and also y multipGer drcuit (MULy) 380. Circuit 375 multiplies the accumulated 
result of the MACx drcuit 370 with the ( l/LE) value to generate the normalized x coordinate of the gamma event ever 
bus 327. Circuit 380 multiplies the accumulated result of the MACy corcuit 365 with the (1/LE) value to generate the nor- 

10 malized y coordinate of the gamma event over bi^ 332. Therefore, the DEP 300 of the present invention computes the 
spatial coordinates (x. y) and the total energy (GE) of each gamma event. Also produced for each event is the peak 
PMT energy (PD) which is the peak signal of the inte^ed channel signals received the DEP. peak PMT address 
(PA) and the focal energy (LE). 

The X. Y. QE and LE values associated with the DEP drcuit 300 of detector 80 are supplied over bus 1220 (see 

IS Rgure 1) to the acquisitfon corrputer 1055. while X. Y, GE and LE values associated with the DEP drcuit 300 of detec- 
tor 80* are supplied over bus 1 222 to the acquisitfon corrputer 1 055. In SPECT mode, the detector pair 80 and 80* per* 
form event detectfon and localization relatively independently because their valid event triggering signals (1240 and 
1 242) do not need to be in ooinctdenca In PET mode, since the event detection needs to be in coincidence, the DEP 
drcuils of both detectors 80 and 80' operate in a form of synchronization t>ased on the coinddenoe of the valid event 

20 trigger signals. 

The spatial coordinates (x. y) of a ganrma event (interaction) are computed by the DEP 300 drcuit using the below 
centroid corrputation: 

Vfa/E^-hWX2*E2+WX3*E3+...4-Wx„*E„ 

25 Ei+E2 + E3+...+ E„ 

Wy/Ei^.Wy2*E2 4-Wy3*E3^....4-Wy„*E„ 
Ei+E2 + E3+...+ E„ 

30 

Where: 

Wxj sthe X weight from drcuit 345 for the i^ PMT of the Ouster 
35 Wy; =4he y weight from circuit 345 for the i^ PMT of the Ouster 
=:the X weight from drcuit 345 for the last PMT of the Cluster 
Wyn y weight from drcuit 345 for the last PMT of the Ouster 
E| :sthe integrated signal for the PMT of the PMT Ouster 
En >^ integrated signal for the last PMT of the PMT Ouster 

40 

The DEP 300 operates as discussed above for each detected gamma event of an imaging session and stores the 
above infbmnation to a computer menxxy storage unit. This informatfon is then transferred to correction electronfos (or 
CPU system) where the data is con-ected for energy linearity, and uniformity This can be done in either the acquisition 
computer 1055 or the image processor 1060. It is appredated that a number of well known methods and drcuitry com- 

45 ponents may be used for collecting the count data output from the DEP 300 drcuHs and for forming an image based 
thereon by correcting the data supplied from the DEP 300 drcuits (e.g.. for nonunifbnnities. ete.) and spatially recording 
the counts. /\ny of these well known methods ntay be used in conjunction within the present invention. 

Figure 3 illustrates the applicat>ility of the atxTve spatial computations and gives an exemplary situatfon. Figure 3 
illustrates a selected PMT dust^ conf iguration as is generated based on the peak PMT address (here it is PMTO) and 

50 based on the PMT address table drcuit 335. The PET imaging node is selected in this example so the PMT duster is 
composed of seven PMTs (six sunounding and one center PMT). The address table corcuit 335 wouM also output a type 
registration to drcuit 340 to indicate that the PMT duster is of a symmetrical or normal type because the peak PMT is 
not located on the edge nor in a comer of the PMT array of the detector head 80. An exemplary x axis 41 0 and y axis 
415 are shown and the weight values for each of the PMTs for the x and y axis are plotted afong the axis for each PMT. 

ss The exenplary event occurs at point 50 within Figure 3 and the arrows extending 0^^ 

li^ received each PMT toward which the associated arrow points. The weight values for each of the PMTs in the x 
and y directions are also dependent on the spatial focatfon of the peak PMTbecause the locatfon of the peak PMT wiU 
alter the PMT type value whfoh is used (in conjunction with the PMT address) to address the weight table drcuit 345. 
Figure 4A and Figure 4B are f tow diagrams atustrating the general processing f taw 460 of the present invention with 
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respect to the operations of a single scintillation detector of the detector pair 80 and 80'. However, the flew 460 is per- 
formed k>/ tx>th detectors. It Is appreciated that these flows pertain to tx)th SPECT and PET imaging modes and that 
valid event trigger detection (block 466) would occur in coincidence for tx>th detectors within PET imaging moda 
Refer to Figure 4A where the procedure enters and receives a signal, over each PMT channel, and after converting 

5 the signal from current to voltage, performs a computer controlled baseline voltage offset for each channel at 462 and 
also performs a computer controlled coarse gain adjustment for each channel at 464. At block 468. each channel is dg- 
itized and also trigger detection is performed at 466 by the event detection logic of the present invention. Trigger detec- 
tion is divided among four area zones across the surface of the each scintillation detector. As discussed above, in PET 
mode, a valid event trigger detectran requires a coinddenoe signal between both d^ectois 80 and 80*. As discussed 

10 above, in SPECT mode, the zonal event detectors utilize a leading edge discriminator where in PET mode the zonal 
event detectors utilize a constant fraction discriminator. Dual integration takes place at blocks 470 and 472 wherein a 
first and a second valid event trigger may be used to integrate separate events occurring close in time. Both integration 
steps output integrated signals to the two stage FIFO circuit and at 474 these data values, per channel, are transferred 
to the DEP sequentially at the completk)n of an integration period. The raw data is sanpled and made accessible to a 

f 5 data processor at 476 and this data is supplied to the calibration table at 478 whteh removes the baseline offset and 
also performs fine gain adjustment of the digital channel signal. 

At 480, the data from the cafibratk>n table is stored for each channel in a buffer. At 482. the present invention deter- 
mines the peak PMT tyy examining the channel data from the calibration table and also at 484 the global energy is deter- 
mined by sunmng the digital data of each channel for the event The peak PMT addre^ 

20 coarse spatial k)cation of the gamma event. At 486, the present inventkxi PMT address table outputs a PMT cluster type 
and also determines the constitution of the PMT cluster based on the peak PMT address {and the selected re80lutx>n, 
e.g., fine or coarse, in one embodiment in SPECT moda) At block 486, smaller sized PMT clusters are used in PET 
mode wNle larger sized PMT clusters are used in SPECT moda 

Referring to Figiffe 4B. the flow 460 of the present inventkm continues at 488 where the circuitry used to peribrm 

25 the centroki computation is reset to initialize for the new computation. At 480. the sequence counter addresses the PMT 
address table so that the first PMT address of the selected PMT cluster is output Rom this value, and also based on 
the PMT cluster type, the present invention generates an x weight QNx) and a y weight (Wy) for the selected PMT 
address. Alsa the buffer 325 contains and supplies the stored integrated signal data for this channel at 494 and at 493 
the dynamk; compressfon circuit outputs the dynamtoally compressed signal value for this channel. At 495. the x and y 

30 multiplicatfon accumulatkxi drculs are used to multiply the weight value times the dynamic compressed signal value 
and accumulate the result for the PMT cluster. At 495. a k)cal energy accumulator also accumulates the local energy of 
the PMT cluster. At 496, the sequence counter inaements and addresses the PMT address tatie for a next PMT 
address until the PMT duster is complete (ag.. the stop indicator of table 335 is reached); flow returns to 490 H the PMT 
cluster is not complete. The above processing then continues with the next PMT address of the selected PMT cluster. 

55 If the PMT cluster is oompleta then flow continues to 497 where tfie x and y multiplicatfon accumulation circuits are 
effectively divided t>y the tocal energy to produce a normalized (x. y) spatial coordinate for the gamma event At 498. 
the pertinent information output from the DEP 300 (including the (x. y) coordinate and the total energy) is output to a 
computer system data processor or a correction boaxd that performs energy, linearity and uniformity correction in known 
manners. When operating in PET mode, a pair of x, y. and energy information is generated from each detector 80 and 

40 80' within a coincidence window. The pair corresponds to two points detected on the scintillatk>n detector pair 80 and 
80' correspondng to two garnnria rays enrvtted in roughly opposite directtons from a po» 

pair is recorded by the acquisition computer system 1055 as corresponding to a discrete PET event and from this pair 
axial and transaxial angles of incidence are computed between the pair. The position of the positron interaction lies 
within the line connecting the two points of the pair. When operating in SPECT mode, the acquisition computer system 
45 1 055 records the detector, 80 or 80', to which the output data for an event applies. 

ACQUISmQN COMPUTER SYSTEM fPIGRAL PROCESSOm 

Refer to F^ure 5 which illustrates components of a general purpose compute system 1112 (whk;h can be imple- 
50 mented within the acquisition computer system 1055) that is capable of executing procedures of the present invention 
for controlling the DEP 300 circuit (ag., control of baseline offset gain, and trigger threshokQ and for performing other 
recited functions. The computer system 1112 corrprises an addressAdata bus 1100 for communk»ting information 
within the system, a central processor 1 101 coupled witti the bus 1100 for executing instructfons and processing infor- 
mation, a random access memory 1 102 coupled with the bus 1 100 for storing information and instructions for the cen- 
55 tral processor 1101. a read only mennoryl 103 coijpled with the bus 1100 for 

for the processor 1101. a data storage devfoe 1104 such as a magnetic or optical cfisk and disk drive coupled wHh the 
bus 1 100 for storing image information and instructions, a display device 1 105 (which can also t>e external such as 
device 1065 of Rgure 1) coupled to the bus 1 100 for displaying information to the computer user, an alphanumeric input 
devtoe 1106 including alphanumeric and function keys coupled to the bus 1100 for communicating information and 
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command selections to the central processor 1 1 0 1 . a cursor control device 1 1 07 coupled to the bus for communicating 
user input information and convnand selections to the central processor 1101 . and a signal generating device ("com- 
munication device") 1108 coupled to the bus 1100 for communicating command selections to the processor 1101. A 
hardcopy device 1 109 (ag.. printer) may also be caused to bus 1 100. 

5 The signal generation device 1 108 Includes a high speed communication port for convnunicating with both DEP 
300 circuits. Input bus 1 120 receives the data output from each DEP 300 such as signals PA 312. PD 317, GE 322, X 
327, Y 332 and LE 337 over buses 1220 and 1222. Bus 1 120 also supplies the raw data output from circuit 310 over 
bus 397 to the processor 1 1 12 (with respect to both DEPs). Output from the processor 1 1 12 are ti^^ 
controlling the fine and coarse baseline offset voltages (e.g.. signals 212 and 214) and the PMT gain signal acljustment 

10 220) for each channel. Processor 1112 also controls the trigger threshold 10& This processor 1112 also controls the 
state of the mode signal 1 252 (Figure 2A3) between SPECT and PET. Processor 1 1 1 2 is also used to program register 
1330 (Rgure 2C) to set the appropriate accumulation interval for both detectors 80 and 80*. The raw data sampled ever 
bus 397 and the control signals generated by processor 1112 may be determined and adjusted in real-time by the 
present invention. 

IS In PET imaging mode, the processor 1 1 12 records together the X, Y, and energy values for a pair of gamma events 
that are detected in coincidence. The pair corresponds to two points detected on ttie scintillation detectors 80 and 80' 
corresponding to two gamma rays emitted in roughly opposite directions from a position-electron interaction. This coin- 
cidence information is used in the image generation procedures that are common to PET scanning. From the coinci- 
dence data of a pair of interactions, a line of intersection of each event can be determined which is projected through a 

20 given plane with the object From these two points, an axial incidence angle aiKi transaxial incidence angle are com- 
puted and recorded. In PET imaging, ottier mathematical manipulations allow a complete mapping of the tfvee-dimen- 
sional distribution of the radionuclide within the organ to be determined. These procedures are weD known in PET 
technology and are not discussed hereia 

The display device 1105 of Figure 5 (or the display unit 1065 of Rgure 1) utilized with the computer system 1112 

2S of the present invention may be a liquid crystal device, cathode ray tube, or other display device suitable for creating 
graphic images and alphanumeric characters recognizatrfe to ttie user. The cursor confrd device 1 107 allows the com- 
puter user to dynanrvcally signal the two dimensional movement of a visible symbol (pointer) on a display screen of the 
display device 1 1 05. Many implementations of ttie cursor control device are known in ttie art including a trackball, finger 
pad, mouse, joystick or special keys on ttie aljphanumeric input device 1 1 06 capable of signaling movement of a given 

30 direction or manner of displacement The keyboard 1 106. ttie cursor control devfoe 1 107, the display 1 105 and hard- 
copy device 1 1 09 comprise the user interface associated witti the image processor 1060. 

DUAL INTEGRATION PER CHANNEL 

35 The present invention provkJes muHvile independent integrators per channel (e.g.. two per PMT) in order to accu- 
rately process high count rates and to effectively reduce problems associated witti pulse pile-up. This feative is partic- 
ularty helpful for processing under PET imaging nxxies. Atthough descrit)ed in a specifk; embodiment utilizing two 
integrators per channel, it is appreciated that ttie present inventfon system may be expanded to encompass a multiple 
tmrbv of integrators per channel (e.g., ttiree, four, five, eto.). As will be (fiscussed in further delaa to fblksw, foe each 

40 integrator added an adcRtionai stage within the serial latch circuit is required. 

Rgure 6A illustrates light intensity response (curve 440) over time for two evente eittter PET or SPECT imaging 
modes). The light intensity response is a well known decaying exponential with a time constant, T. Event 1 occurs and 
decays over time as shown by curve 440. Due to ttie characteristics of the crystal 81 used in ttie detector 80, at ttie end 
of ST time intervals, most of ttie usable light intensity is radiated. However, ttiere is some time value R, less than 5T, ttiat 

45 can be used and gives a suffk;ient amount of energy to register event 1 (e.g.. energy associated witti region 430). Prior 
art systems will attempt to utifize this deaeased amount of en^gy (regfon 430) for registering evente during periods of 
high count rate (e.g. , if a second event occurs before ST as shown as cun^e 435). However, ttiis is not advantageous 
because as ttie separation between evente decreases below ST, less and less event energy is captured and spatial 
compulations become less accurate. Furttier, at some point. (e.g., less then R) ttie temporal separation b^een two 

so evente wOl t)ecome too small and neittier event can be registered. At higher count rates ttie spatial accuracy of the prior 
art systems deaeases signif k»ntiy. 

The present invention, on the ottier hand, provides a mechanism for integrating the energy over 5T (for example) 
for botti event 1 and event2 because dual integrators are supplied per channel for all the PMTs of the array. Therefore, 
one integrator may sample and integrate the light intensity for event 1 over response 440 for 5T duration (or any ottier 

ss programmable duration) and ttie ottier can sanple and integrate ttie light intensity for evente over response 435. During 
periods of hig^ count rate, ttie present invention does not sacrifice energy intensity when sampling each event when 
integrating over two evente that closely occur in time. Furttienmore, since two separate integrators are used, the present 
invention can accurately register two evente even though ttieir temporal separation is less than R. The present invention 
tiierefore provUes higher accuracy at higher count rates over ttie prior art 
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The present invention dual integrator embocfiment utilizing two independent integrators is now discussed. As dis- 
cussed with reference to Rgure 20. the integration circuitry of drcuit 280(i) for a given channel includes two integrators 
238 and 240 each having independent accuniulators and each coupled (via a mux 241) to a two stage sample and hold 
drcuit Oatch 242 and latch 244). Each integrator is independently triggerak)le and separately and independently inte- 

5 grates its assodated channel signal over a programmable period. Trigger signals are transmitted over kxjs 130 and 
when received, act to reset the accumulator of an idle integrator. At the end of the progranvnaUe integration interval 
maintained by regster 1 330 (Rgure 2C). tor an event the present invention transfers the data of the second stage (244) 
to the DEP 300, moves the data of the first stage 242 to the second stage and moves the data of the accumulator of the 
completed integrator into the first stage of the hold drcuit This way. both integrators may be sampling, simultaneously. 

10 different events. They can each be independently triggered and at the end of the programmable sample period, the 
accumulator stores its result in the two stage hold drcuit 

The akxTve procedure operates most accurately when the two integrated events are suffidently separate from one 
another such that their energy dissipation across the detector does not overlap or interfere in the integration computa- 
tion. Fbr instance, refer to Rgure 7 which Dlusbates an exemplary PMT array. A first event occurs over PMT 7. therefore 

IS the PMT duster 75 is corrposed of PMTs 1 . 8, 20. 1 9. 36. 1 8. and 7. The trigger pulse resets integrator 238 (per chan- 
neO which then integrates the energy for the first event fbr all 55 channels. Before the integration is complete for this first 
event a second event occurs over PMT 38. so the PMT duster 73 is composed of PMTs 38. 37. 22. 23. and 39. Inte- 
grator 240 (per channel) is reset and integrates the energy associated with the second event fbr all 55 channels. Since 
PMT duster 75 and 73 are suff idenfly separate, the contribution of energy associated with the second event over chan- 

20 nels 1. 8. 20. 19. 36, 18 and 7 is relatively small and does not interfere with the integration computation fbr the Rrst 
event Likewise, the oontrOxjtion of energy assodated with the f irst event over channels 38. 37, 22, 23, and 39 is rela- 
tively smalt and does not interfere with the integration ootrputation for the second event. 

At the end of ttie computation for the first event the integrated channel signal data of latch 244 is output to the DEP 
300. fbr a tfven detector, the infected channel dgnal of latch 242 is output to lat^ 

2S later of inte^tor 238 is output to latch 242. At the end of the computation for tiie second event, the integrated channel 
signal of latch 244 is output to tiie DEP. the value of latch 242 is output to latch 244 and the value of the accumulator of 
integrator 240 is output to latch 242. The first and second events will be processed t>y the DEP 300. 

Figure 6B illustrates in more detail the process perfamed t>y the present invention to perform dual integration per 
channel . Each channel of each detector performs this process. The process shown 501 represents the process fbr inte- 

30 grator A but it is appreciated that the process for B (block 516) is identical except tiiat the "B" control signals are used. 
As shown, the process starts at 510 in response to a valid event trigger signal, Start(t0). Then, ttie A accumulator is 
deared (ag., integrator 238) by assertion of a CLfWX^A control signal and at 512 this accunulator is enabled to inte- 
grate by assertion of an ENACCA control signal. At block 514. if Start(t1} is detected ttien at tkxk 516 ttie process for 
integrator B is started and operates concurrentiy witti process 50 1 . If not, then at bk)A 

3s maUe duration integration fbr A is completa If so not ttien the process returns to 51 2 where integration continues. 

When integration for A is complete over the programmable interval, at btock 520 a control signal fbr MUX 241 
selects ttie data from integrator 238. At block 522 FIF01 (242) is checked if empty and if empty, at bk)ck 524 ttie data 
Is latched into FIF0 1 . At block 526 if RF02 (244) is empty ttien data from FtFOI is latched into FIF02 at block 528. At 
block 530 ttie data transfer from FIF02 to the Di^ Event Processor is started fbr aD channels. During ttie period from 

40 block 51 0 to block 522 ttie A integrator is busy. It is appreciated ttiat ttiere is only one valid event trigger signal and it is 
classified as 8tart(t0) or Start(t1 ) by its tenrporal relationship to ttie ottier trigger signal. During ttie A busy period, a trig- 
ger signal will cause ttie B process to start The regster 1330 (Rgure 2C) can be loaded with the programmable inte- 
gration period at any time but is typically loaded before an imaging sesskxi is started and is not altered during ttie 
session. 

45 

VARIABLE PMT CLUSTER CONSTITUTION 

As discussed above, the PMT address table drcuit 335 of the present invention contains a lookup table that pro- 
vides the addresses of the group of PMTs ttiat constitute the PMT duster fbr a given event based on the peak PMT for 

so that event (supplied from drcuit 320) and based on a count value sillied over bus 357 from ttie sequence counter 390. 
A centroid is computed (using a centroiding computation) based on ttie PMT duster to arrive at a spatial coordinate 
value of the event. In this way. the PMT duster constitution of the present invention varies for each peak PMT. According 
to ttie present invention, also associated with each PMT of the detector array is a type registration or classification that 
desabes the type of PMT duster that is formed. 

ss Due to the Ngh count rate experienced during PET imaging, ttie present invention, via ttie computer system 1112. 
generally controls ttie size of a PMT duster to be no more ttian seven PMTs per cluster while in SPECT mode ttie typ- 
ical PMT duster size is allowed to indude more PMTs, such as 19 PMTs per duster. The exact nunrt)er of PMTs per 
duster is not essential, however, ttie present invention advantageously allows duster sizing to be reduced when oper- 
ating in PET mode while PMT duster sizing is increased during SPECT mode. Therefore, during initialization of a PET 
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imaging mode of operating, the system 1112 automatically controls the programming of the present invention to select 
and utilize a reduced number of PMTs per duster while during Initialization of a SPEC! imaging mode of operating, the 
system 1112 automatically controls the progranvning of the present invention to select and utilize an enlarged number 
of PMTs per cluster. Spedfically. to initialize under either PET on SPECT imaging nKxJe, computer system 1112 loads 
5 a predefined PMT duster table (Figure 8) into circuit 335 so that the proper PMT cluster sizes are availabia The prede- 
fined tables (one for PET and one for SPECT) can be stored inunrt1102or1103or1104, for instance, before down- 
loading Into drcuit 335. 

With regaid to variable duster shape, the present invention provides four different types of exemplary PMT dusters: 
1 ) normal; 2) long edge; 3) short edge; and 4) comer. The normal PMT duster anticipates the peak PMT to be located 

10 In an area of the PMT array that can bG surrounded by other PMTs such ttiat the PMT duster is substantially symmet- 
rical about botti axis. Such a PMT duster is shown in Rgure 7 as PMT duster 75 wHh PMT 7 as the peak PMT which 
is surrounded k>y outer PMTs. The determinatx>n of the peak PMT gives a coarse spatial location for ttie event and the 
present invention PMT address table 335 is able to utilize this coarse location of the event in determining the PMT dus- 
ter constitution for that evem so thai the fine spatial location niay be provided vi^ DEP 

IS 300. 

Two PMT duster types (long edge and short edge) correspond to PMT dusters having peak PMTs kx;ated on tiie 
edge of the detector PMT array. Edge PMT dusters are not synfmetrical about one axis. The PMT array of the present 
invention is rectangular and contains a kxig edge and a short edge and depending on the edge location of the peak 
PMT, the PMT cluster defined thereby can be a long edge PMT duster or a short edge PMT duster. The PMT dusters 

20 dl these types are different because the light distribution at the tong and short edges are cfiflerent An edge type PMT 
duster is shown in Figure 7 as PMT duster 73 having PMT 38 as the peak PMT The fourth type of PMT duster of the 
present invention is a corner PMT duster and the PMTof this type is k>cated in a comer type of the PMT array, such as 
PMT duster 71 of Figure 7 having peak PMT 46. This PMT duster type is not symmetrical about either axis. It is appre- 
ciated that the above PMT cluster types are exemplary and many other PMT duster types may be utilized within scope 

2s of the presem invention depending on the particular geonrietry of the PMT array utilized £^ 
tor head. 

It is appreciated that for a given peak PMT address, the type value for the PMT duster can change depending on 
the selected resolution of the spatial computation. For instance, a normal type PMT duster at high resolution having the 
sanrie peak PMT nriay be different from the PMT cluster generated al low resol^^ 

30 be an edge type duster). This is the case because at least in SPECT mode, at higher resolution, more PMTs are 
required (e.g.. 17 or 19) to complete the PMT duster rather than 7 fbr the low resolution duster and these additional 
PMTs may overrun the edge of PMT anay. 

Type fields are important because the PMT duster type will effect the x and y weights assigned to a particular PMT 
of a particular PMT duster in the spatial oornpulalion. R)r instance, PMT cluster types that are symmetrical about only 

3s one axis (ag. , edge types) will have modified weight values assigned to those PMTs of the PMT duster that are k)cated 
along tfie axis ttiat is not symmetrk»l. For instance, refer to Rgure 3 whteh illustrates a spatial computation based on a 
normal PMT duster type. The spatial coordinate Is corrputed from an average of a sum based on the weight of a PMT 
multpGed by the Integrated channel signal of the PMT for each PMT of the PMT duster. The weig^ values fbr the PMTs 
must be adjusted in the computation of a spatial coordinate ak)ng an axis for which a PMT duster is not syirwnetrical. 

40 Referring to Figure 3, assume that PMTs 2 and 3 were not available because PMT 0 was k)cated along an edge. The 
resulting PMT duster is composed of PMTs 0. 1 . 6, 5, and 4 and is not symmetrical about the X axis. For ttie coonfinate 
computation of the X axis coordinate, the average cak;ulation woukJ be skewed or shifted to the left t>ecause the PMTs 
of the right (e.g. . PMTs 2 and 3) are missing. Therefore, the weight values of the PMTs of the resulting PMT duster must 
be reduced to compensate for the skew to the left. The same is true for corner PMT dusters* however the weight values 

45 must be adjusted fbr the compulation of both coordinates because comer PMT dusters are not symmetrnal to the X or 
Yaxis. 

The weight adjustment of the PMTs based on the PMT duster type f iekl is made by the weight drcuit 345 and will 
be discussed further below. 

Refer to Rgure 8 which illustrates the format of the PMT address table drcuit 335 fbr low resdution selection In 
so SPECT nrKxIe or fbr PET nrxxie. The PMT address table drcuit 335 is addressed by the peak PMT address value and 
also addressed by the coum value (here shown from zero to n). Circuit 335 conteins 

of tfie detector array. For each peak PMT, the drcuit 335 outputs a PMT cluster type value and the PMT addresses of 
the PMT di^er. The output of the drcuit 335 is (1 ) a PMT duster type value and (2) the PMT addresses (a 'PMT MsT) 
of the PMTs of the PMT duster defined by the peak PMT address. Since the PMT dusters are programmable and of 
ss variable size and PMT number, a "stop" incficator is placed at the end of the PMT list (or included as part of the last PMT 
address entr^. Exemplary date is shown in Rgure 8 and corresponds to the PMT dusters shown in R^^e 7. The first 
entry shown of Rgure 8 relates to PMT duster 74 of Rgure 7 and PMT 7 is the peak PMT and the PMT duster is a 
normal type having PMTs 1 ,8.20,19.36, and 18. Entry 38 of Rgure 8 relates to PMT cluster 73 of Rgure 7 and is an 
edge type PMT cluster. Entry 46 of Figure 8 relates to PMT duster 71 of Figure 7 and is a corner PMT duster. 
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Ihe data stored in the memory cfrcuit 335 that is used to formulate the PMT clusters for each peak PMT can be 
programmabie and may t>e downloaded from the conrputer system 111 2 at initialization of an imaging nrxxJe or between 
SPECT/PET selections. In such an embodiment, different datasets may be loaded into the memory circuit 335 for dif- 
ferent configurations. Altematively the circuit 335 may be implemented using static ROM memory. In one embodiment, 
5 two separate PMT cluster tables (Rgure 8), one for PET and one for SPECT imaging, can be loaded into circuit 335 
with a mode selection (line 1252) being fed into drcuit 335 for selection between the two tables for PET and SPECT 
moda 

Based on a dock signal, the sequence counter 390 of the present invention will present the count i'lM ever bus 367 
(one at a time) and this count value will address the drcuit 335 along with the peak PMT address to output (1) the PMT 

10 duster type and (2) the PMT address of the PMT duster as shown in Rgure 8. 

It is appreciated that depending on the desired resolution (or count rate) of the gamma camera, or the operational 
mode of the camera (e.g, PET or SPECT), the PMT address table 335 of the present invention will output different sized 
PMT dusters for each PMT duster typa For instance, if high resolution spatial determination is required, or if SPECT 
imaging is being performed, or low count rate is expeded, then the PMT dusters wfll be increased in size to include 1 7 

IS to 1 9 PMTs for a normal PMT duster Cmstead of seven for normal duster types in lew resolution). As discussed above, 
if PET imaging is desired then 7 PMTs per duster are supplied to allow high count rate with minimum pileupi Edge and 
comer PMT dusters will be increased accordingly in number. Therefore, in an alternative embodimerrt of the present 
invention, the PMT address table 335 receives an additional signal indicating fow or high resolution and this signal will 
address the table to supply the appropriate centroiding information based on the required resdution. Altematively. the 

20 entire table 335 may be retoaded with a difl^em data set to vary the resolution in SPEC^ mod& 

Given the design of the PMT address taUe of Rgure 8. the number of PMTs of a given PMT duster may easily be 
increased since the stop indicator, which marks the completion of the PMT duster, is readily acQusted. Further, as 
stated above, the PMT duster type corresponding to a particular peak PMT address may vary from low to high resdu- 
tion settings, or fbr PET or SPECT imaging modes. 

2S 

VARIABLE PMT WElGhfTS PER PMT 

According to the present invention, the weight table drcuit 345 oulputs x and y coordinate weigfits for each PMT of 
the PMT duster based on the PMT address arvj tfie PMT duster type information, both of whk^h are generated by the 

30 PMT address table 335. Rgure 9 illustrates the weight table drcuit 345 of the present inventfon. Depending on the type 
of PMT cluster (e.g., normal type, long edge type, short edge type , comer type, or other type) that the PMT is contained 
within, the x and y coordinate weights output from the drcuit 345 for the PMT will vary. For each PMT address (ag., 
from PMTO to PMT54). the present invention provides a separate and programmable weight value for each coordinate 
compulation (ag., Wx and Wy) that varies PMT duster type. The values Wx are output over bus 362 and the values 

35 Wy are output over bus 357. As will be (fiscussedbdcw, the determinalion^ 

tion for the evem and the present invention wdght table 345 is able to utilize this coarse focatkm deter- 
mining the proper weight values to assign the PMTs of the PMT dieter. The fine spatial focation is computed via the 
centroid computation of the DEP. In such manner, the PMT address value signal over bus 372 and the type signal over 
bus 377 address the drcuit 345. 

40 Since some PMT dusters are not synrunetrical about a given axis, for instance ft^ 

PMT cluster, the weight values associated with these axis are adjusted or varied in order to balance out the coordinate 
computation. This is accomplished t>y the present invention in order to compensate for the missing PMTs (of one axis) 
that are not available to provkie a symmetnc computatioa For comer PMT dusters, the weight values associated with 
tx>th axis are adjusted to compensate for the missing PMTs (of both axis) needed to provide a symmetrk; computation. 

45 Typically the weight values are decreased for certain PMTs in order to perform the above balandng. The values weight 
labfo 345 of the present invention may also be devdoped empirk»lly bas^ 

Further, other factors such as the crystal boundaries, optical interfaces, and PMT photocattKXle properties can 
make the PMT contrfoution different depending on ttie location of tiie event Since, the peak PMT address gives some 
indication of the coarse focation of the event, the weight taUe 345 can conpensate for the abm 

so variable wdghts. 

Therefore, the present invention provkJes the ability to adjust or alter the weight values for a given PMT depending 
on the PMT duster type in which tfie PMT is located. Depend ng on the location of the peak PMT, the weight values for 
the PMTs used in the centrokJ computation may alter. These weight values are also deperxient on tfie peak PMT 
address since the peak PMT address defines the PMT duster type within the present invention. The atNlity to assign 
55 different weighting factors based on the peak PMT location permits higher accuracy in the centroid computations and 
reduces the demands on the oonrection processing steps. This contributes to alfowing the detector 80 to have larger 
field of view without increasing the crystal dimensfons. 

The data stored in ttie memory drcuit 345 of Figure 9 ttiat is used to provkie the variable weights for each PMT can 
be programmable and may be downfoaded from the computer system 1 1 12. In such an embodiment, diff^ent datas^ 
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may be loaded nnto the meifiory circuit 345 for different oonffigurations. Altematively, the circuit 345 may be implemented 
using static ROM memory. 

in operation, as the counter 390 counts sequentially, the PMT address table 335 outputs a sequential listing off PMT 
addresses within, the PMT cluster. The type signal generated by Ixjs 377 addresses the MSBs of the of the memory 
5 345 and the PMT addresses are the LSBs. As each PMT is generated over bus 372, the memory circuit 345 generates 
an X and Y weight value (over buses 357 and 362) associated with the current PMT output over bus 372. This informa- 
tion is fed to the centroicfing circuitry for computation of the spatial coo^ 

AUTOGAIN CORRECTION 

10 

As will t>e discussed further below, the autogaln calS>ration of the present invention contain two phases: (1 ) an initial 
calibration; and (2) a routine calbration. The routine calibration can be performed with the colfimator installed or with 
the collimator removed. Autogain correction, to the extent it operates in conjunction with a collimator is applicable only 
to SPECT imaging modes of operation because the PET imaging modes do not require oollimator& 

IS The present invention provides a procedure for automatically calikrating the gain of each PMT channel without 
removing the collimator during the routine gain calit)ration phase. This is benefictal at least since the gain calibration 
may be performing during periods of quality control wrttKxit removal of the collimator while imaging a sheet source (ag., 
Co-57) and thus saving time and effort Further, since collimator removal is a laborious and time consuming procedure, 
the automatic gan caGbration of the present Invention provides all efficient mechanism for gain connection. The oompu- 

20 ter system 1112 adjusts a preamplif ication gain (C^) for each PMT by using coarse and fine gain values. TTie coarse 
gain adjustment is accomplished tyy generating a digital gain value that is converted to an analog signal by a DAC and 
applied at input 220 of the preairplif ication circuit (280(i)) for each channe](i). The fine gain adjustment is implemented 
in the calibration tMe 315 as a looK*up tabia The computer system 1112 stores tiie current value of both the coarse 
and fine gain (G^) for each channel. The computer system 1112 also perfbms the automatic gain adjustment proce- 

25 dure as will be discussed in further detail beHcw. 

The effective gain associated with each PMT channel has two separate components. The first component is the 
physical gain associated with the individual PMT itseff (Qt) and this gain is established by the physical characteristics 
of the PMT. The second component is the corresponding preamplifier gain {Qp). As is known, the PMT pfiysical gain, 
Gt, may vary over time and change in the long term (e.g., over hours or days). These variations in the Qt create Itie 

30 problem fliat the present invention solves by adjusting ttie Qp gain associated with each PMT channel individually to 
compensate for the changes in Gt toobtain a stable (fixed effective gain) for each PMT cfiannel. 

The preamplif ication gain, Gpy contains two components. For each channel, the first component of the preamplifi- 
cation gain Gp is the gain applied to the output cfiannel si^l within the preamplif ication stages, 280(i), for each chan- 
nel (e.g., coarse adjustment); tiiis component is applied at gain circuit 222 via input 220. A second component (ag... 

3S fine adjustment) of the preampliffication gain Qp is supplied by the calibration table 31 5. if the gam of a particular chan- 
nel does not vary by more than some small percentage (e.g., 5%) then the gain value found wrihin the calbration table 
for that channel is altered, not the gain applied at the preampliffication stage for tfiat cfiannel. 

The efffective gain, QtOp, is to remain f ixed over time to yield the same output signal off a particular racfionudide 
emitting pfKiton with distinct energy. The fixed gain of each PMT cfiannel will give stat)le x, y ooonfinates and total 

40 energy values which will increase the time t>etween required recalit)ration of tfie camera's energy, linearity and uniform- 
ity correction lectors. Another advantage of the autogain caEkXBtion of tfie present invention is that it t)rings bodk, tfie 
camera to its state when it was initially calit)rated and when uniformity correction factors were generated. 

In view of this, since the gain Gt may vary due to the physical cfiaracteristics of the PMT, the present invention com- 
pensates for this variation by altering the computer controlled value of the gain that is applied to data from each 

46 channel. In this way, the automatic gain correction procedure of the present invention will maintain the resultant gain 
(ag.. Gt * Qp) at a fixed value for each channel. The actual value of this fixed amount can be detemwied and/br meas^ 
ured during manufacturing or at the site of Installation and operation of the scintillation detector. Since energy, linearity 
and uniformity corrections are based on a fixed gain of each of tiie calforated photomuttiplier cfiannels of the array, it is 
advantageous to maintain this calibration throughout the operational cycles of tfie gamma camera system so tfiat the 

50 gains of tiie PMT channels are dosely calibrated to match the gain when the energy, linearity and uniformity correction 
factors were calibrated. 

Central and Ot)scured PMTs. Rgure 10 illustrates an exemplary PMT array of the present invention as well as a 
collimator 83 that is positioned in front of the PMT anray. The collimator 83 contains an inner region 61 5 composed off a 
honeycomb of holes witti lead septas allowing racfiation of a particular incident angle to pass there through. The cdli- 
55 mator 83 also contains a second region (edge region) 83 that is solid lead and partially and fully covers some edge and 
comer PMTs. PMTs having their surfaces located totally under region 615 are the central PMTs (centrally located) and 
PMTs having a portion of their surfaces located under region 61 0 are called obscured PMTa The present invention pro- 
vides an automatic procedure and apparatus for calibrating tiie gain factors of each PMT of the array Cinduding those 
PMTs that are obscurecQ while leaving tiie collimator installed during the routine cdibration procedure. Since PET 
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imaging does not utifize a collimator, PET imaging does not utilized obscured PMTs. 

Each PMT has an area 570 associated with tt directly ever the PMT surface, see region 570 of PMT 24 for 
instanceOn Rgure 10). Although only one area 570 is shown, it is appreciated that each PMT of Figure 1 0 has an asso- 
ciated region 570 located above Its individual suriace The actu^ 

5 it may circular, square, hexagonal, eta). In the preferred embodiment of the present invention, this area 570 is drcUlar 
and covers substantially all of the area above a given PMT. TNs central region 570 will be used for calibration of the 
preamplification gains (including adjustment of the calibration table) of each PMT when the collimator is removed and 
each centrally located PMT when the collimator is installed. (Some PMTs can be located only partially over the Nal 
crystal layer and therefore will receive less light) 

10 Strip Regions 910. Rgure 1 1 illi^tes specialized regions 91 0 (one associated with each obscured PMT) utilized 
by the present invention for performing gain calibration for the obscured PMTs of the scintillation detector. Rgure 1 1 
illustrates two obscured Pf^s 39 and 40 that are obscured by the collimator's solid lead region 61 0 (when installed) as 
well as a group of centrally located PMTs (0-6) that are situated inside region 615 where the collimator consists of holes 
pemnitting radiation to pass through. The region 91 0. or area is associated the obscured PMT 39 but a separate 

IS strip area 91 0 is also assodated with each obscured PMT of the scintillation detector. Rgure 1 1 illustrates only one 
such strip re^on 91 0 assodated with obscured PMT 39. This region 91 0 may be qif a variety of different geometries and 
in the preferred embodiment is semi^nular and extends in an arc fashion sunrouncfing a region that is located substan- 
tially equidistant from the center of its corresponding obscured PMT 39. 

It is appredated that this strip region 91 0 may vary in size and shape consistent with the present inventioa How- 

20 ever, in a preferred embodiment of the present invention, this region 910 is located such that it extends substantially 
equidistant from the center of its assodated PMT and generally adopts a thin shape as shown in Rgure 1 1 . A further 
rec^irement of strip 910 is that it extend substantially (if not fully) within the region 615 of the ooOimator that allows 
gamma rays to pass ttirough. 

When the cdRmator is InstaHed onto the sdntillation detector, the present Invention is not capable of perfbrnning 

25 gain calibration on the obscured PMTs utilizing an area located dtrectty over their surfaces, such as region 570. 
because no gamma events are detected under region 610 due to the presence of the collimator's edge. Therefore, ttie 
present invention perfonns gain calit)ration for obscured PMTs by recorcfing information assodated with gamma events 
that are detected within the exposed strip area 910 assodated with each obscured PMT. In general, the present inven- 
tion measures the integrated channel response of an associated PMT. e.g.. PMT 39. in response to ail gamma events 

30 that occur wHhin the strip area 910 over a sample intend during the initial calibration. This information will be used to 
adjust the gain of ok>scured PMT 39. 

PMTs of a PMT duster detect different amounts of energy responsive to a gamma interaction depending on their 
spatial relationship to the interaction. Refer to Rgure 12 which illustrates an exemplary and typical PMT cluster 553. R>r 
a g^en gamma interaction (ag., 140 keV), approximately 4000 separate 3eV visft)le sdntillation photons are emitted 

35 from the crystal 81 and of which only 50% are detected. Depending upon the size of the PMTs^ if the gamma interaction 
occurs within the center of region 570. the center PMT will coOect and report approximately 45 percent of the detected 
light energy. The surrounding adjacent six PMTs each cdlect approximately 8% of the detected energy. For a nurrber 
of detected gamma events located near the center of region 570 of ttie center PMT. the center PMT should produce an 
integrated output voltage proportional to approximately 1000 3eV light photons. Each of ttie acijacenfly sunrounding 

40 PMTs shouM output a Integrated nraan voltage proportional to 160 3eVligmpho 
within 570. 

Refer to Rgure 13A, Rgure 13B. and Figure 13C which illustrate processing steps executed by the computer sys- 
tem 1 1 1 2 of ttie present invention for perfbrnning automatic gain calbration/^justment of ttie prearnplif ication gains Qp 
(ag.. assodated witti drcuils 280(i) and ttie adjustment wittiin ttie calibration table 315) for each channel of the sdntil- 
45 lation detector. 

Figure 13A IHustrates ttie process flow 950 of ttie presem invention ttiat is used to perform inW^ 
tion of the sdntillation detector and is typically perfomied prior to ttie energy, linearity, and uniformity calibration. The 
energy, linearity, and uniformity calixation can be performed at the manufacturing site or can t>e performed at the oper- 
ational site of the camera system. This is done prior to operational imaging of a patient As shown, at step 957, the cd- 

so limator is removed (if present) from ttie sdntillation detector and ttie array of photomultipliers is irradiated witti a unifbrm 
flood field of known energy gamma radiation (ag.. a known isotope Is utilized). I=br all events ttiat are detected wittiin 
the regions 570 of Rgure 1 2 (for each PMT). ttie total energy of the events are recorded, individually, irito memory 1 102 
and assodated witti ttie appropriate PMT such that a distribution is formed for each PMT of only ttiose events that were 
detected above ttie PMT. The memory 1 102 stores the information in a matrix form that associates a given PMT witti a 

S5 distribution of ttie total energy of gamma events detected wittiin regton 570 fbr that given PMT. At block 957. a meas- 
ured peak total energy value of ttie distribution is determined by ttie computer system fbr eaich PMT 

Since a known isotope is utilized, its total energy is known (e.g.. nominal total energy peak). At block 959 for each 
PMT. ttie measured total energy peak value, as reported by ttie sdntillation detector at block 957, is compared by ttie 
present invention against ttie known ("nonrinaT) total energy peak value assodated witti ttie isotope. Wittiin block 961 
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of Figure 13A. the gain of each PMT channel is then adjusted such that the measured total energy peak value matches 
the known or nominal total energy peak value. 

For instance, since the energy detected by the scintillation detecta responsive to each gamma event of the flood 
f ieU is kncwn (e.g., X), the peak energy catoulated and associated with all PMTs is compared against the known value 

5 (X), and the gain (G^) of the PMT is adjusted up or down accordingly such that the measured value and the nominal 
value are equal. It is appreciated that both the gain values associated with the preamplification circuits 280(1) and also 
associated with the calibration table 315 are ac$usted at block 961 . The gain values associated with the preamplif ication 
circuits are a ooaise gain and the gain values of the calixation table 3 1 5 are fine adjustment values. At bkx:k 961 . the 
gain values are also appfied to the preampGfication channels and the calibratk)n table 315 for each PMT channel. 

10 Refer to Figure 12 which illustrates a PMT cluster of PMTs 0. 2, 10, 1 1 , 12, 4, and 3 (center). The total energy of 
events occumng within region 570 of PMT 3 is recorded and stored over a numt)er of events and a measured peak 
energy value is determined from the energy distnbution (at block 957). For events occurring within region 570 for PMT 
3. PMT 3 is inherently the central PMT (of the PMT cluster). Therefore, it is expected that the central PMT detects 
approximately 45% of the li^ energy ()0 detected tyy the detector. If the total energy of the detected event (at region 

IS 570) is 5% higher than the expected value of X. then the gain of the center PMT 3 is adjusted downward by 10% (at 
block 961) because approximately 45% of 1 0% adjustment will equally compensate for the detected 5% deviation. This 
is performed for each of the 55 PMTs of the PMT array. 

However, the above gain (Op) correction ignores the contribution of the peripheral PMTs. It is true that one or more 
of the peripheral PMTs (ag., 0, 2. 1 0, 11 . 1 2, or 4) may have a deviant gain which, for example, may cause a 5% energy 

20 deviation for events detected within 570 of PMT 3. If this is the case, then the above correction to the gain (G^) of PMT 
3 may have been made in error. Therefbre, the present invention, at block 963. performs the above gain calibra- 
tionAadjustment for each of the 55 PMTs a number of different times in an iterative procedure (ag., 1 0 to 20 times but 
the nunber is programmable) to account for the possibility of inaccurate adjustments. Based on this iterative procedure, 
the effects of inaocirate gam acljustments become signif icanfly reduced as each PMT is eventually adjusted as a center 

2S PMT lor each iteration. At the end of the processing of the last iteration, the final preamplification gain (G^) assigned to 
each PMT is recorded in memory 1 1 02 (or other storage unit) at step 965. Each gain value is stored associated with its 
PMT 

Referring to Figure 1 3A, block 967 is then entered. The collimator remains removed from the detector and the PMT 
anray is again Qluminated with the same flood field radiation. For each "obscured" PMT (ag.. not obscured during block 

30 960, but obscured when the coDimator is subsequentiy installed), a region 910 is defined that extends into the field of 
view of the detector that will encompass the open region 615 of the collimator (when installed). For each 'obscured' 
PMT, each gamma interaction that occurs wHfiin its region 910 Is recorded by the computer system 1 1 12 at ttock 967. 
The integrated cfiarmel output for tfiat "obscured" PMT is then recorded to memory 1 1 02 and an energy distrilxjtion is 
formed wherein the peak energy value Is then determined at block 969. In the preferred embodiment tfie "obscured' 

3S PMT is a peripheral PMT by definition for all gamma events that occur within its associated region 91 0, therefora the 
"obscured' PM is expected to receive a fraction, C, of the total event energy (wherein C is less than 10 percent). For the 
exanple of Rgure 11, tfiis would represent a fraction, C, of the value X (2000 photons of 3eV each) as reported kiy the 
detector. The measured peak integrated channel signal detected at block 960 for the "obscured" PMT should register 
roughly a fraction of X. or C"X. 

40 For instance, (see Figure 1 1) for all events occurring within the region 910 associated wifli PMT 39, the integrated 
channel signal ou^xit associated with PMT 39 Is recorded until a cfistribution is formed In memory 1 1 02. This is per- 
formed for each ot>scured PMT (with respect to its associated strip area 910) over a number of gamma interactions 
(ag., 500, but the number is programmable). The peak PMT signal of the distribution of integrated channel signals for 
each obscured PMT is then recorded into rnenmy 1102 or other storage unit associated with ^ Processing of 

45 block 950 then retuma 

At the completion of the process 950, the present invention calbrates each PMT tube of the array according to 
events that occur over each PMT and records this gain value in memory Process 950 also records the peak integrated 
cfiannel signal for each 0k)scured PMT in response to events that are detected wHh the ot)scured PMTs Individual strip 
region 910. 

so Refer to Figure 13B which illustrates the routine calibration process 970 of the automatic gain correction of the 
present invention when the collimator is removed. Routine calibration 970 can be peribrmed at the operational site for 
calft)ration of the sdrrtillation detector while the collimator is removed. The processing of blocks 971 -979 are analogous 
to the processing tasks 957-965 of Figure 13A. Therefora at block 979 of Figure 13B, the corrected gain values are 
stored in memory 1 1 02 and implemented within the preamplification stages 280(0 and the calibration table 315 for each 

55 channel. Since the coDimator is removed, there are not any actual obscured PMTs. 

Refer to Figure 13C which illustrates the routine calibration process 980 of the automatic gain oon-ection of the 
present Invention when the collimator is installed. Routine calixation 980 can t>e performed at the operational site for 
calbration of the sdntillation detector while the collimator installed. At block 981 . a uniform flood field source is placed 
in front of the oollmfiator. Interactions occur over the PMT array and the total energy of events detected cv^ 
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region (570) of each central PMT (ag., non-cbscured PMT) are recoided in memory 1 102 and associated with that 
PMT. A distrftxjtion is recorded for each central PMT and a measured total energy peak value is determined for each 
central PMT. Block 981 is analogous to block 957 except the processing for block 981 is applicable only to central PMTs. 
The presence of the collimator does not interfere with this process because the centrally located PMTs are under the 

5 open section 61 5 of the collimator. 

At btock 983, for each obscured PMT the present invention records the integrated charviel signals from the 
obscured PMT only for gamma events that occur within that ot)scured PMTs associated strip region 910. A separate 
cfistribution of Integrated channel signals is recorded for each obscured PMT and stored in memory 1 1 02 (or other stor* 
age device). As slated above, the strip regions 910 associated with the obscured PMTs are aligned with the central por- 

10 tion of the collimator. The processing of bfock 983 is analogous to block 967, except at block 983 the collimator remains 
installed. It is appreciated that bfocks 981 and 983, of process 980, can occur coinddentally 

At block 985 of Figure 13C. the measured total energy peak values of the central (e.g., non-obscured) PMTs as 
reported from block 981 are compared against the known (nominal) total energy peak values for each central PMT 
channel. The gain value for each central PMT is then adjusted acoordin^y. The processing off btock 985 is anatogous 

IS to the processing of bkx:k 959 except only central PMTs are processed in block 985. The new oom(^^ 
the central PMTs are sk)red in menrwry 1 102 (or similar storage device). 

At block 987 the present invention, for each ot>scured PMT channel, compares the measured spectrum peak value 
(that was measured within block 983) against the stored peak value associated with that obscured PMT channel (that 
was stored by btock 969). A! btock 987, the present invention computes a new gain (tsy increasing or decreasing the 

20 current gain) for each obscured PMT such that the measured value matches the stored value fbr each ol>scured PMT. 
This can t>e performed by a straight ratfo computation (ag.. if the measured spectrum peak value reported from block 
983 is 1 0% larger than the stored value from block 969, then the current gain for tfiat obscured PMT is reduced by 10%, 
eta). At block 987, the newly computed gain fbr each obscured PMT is then recorded in memory 1 1 02 (or similar stor- 
age device). 

25 At btock 989 of Figure 1 3C. the present inventton then applies the newly computed gains of: (1) the central PMTs 
as computed in btock 985; and (2) the obscured PMTs as computed in block 987, to the appropriate preamplificatton 
circuits 280(i) and the appropriate locations within the calit)ration table 315 for each channel. The preamplification dr- 
cuit contains the coarse gain control while the calibration table 315 prcvkie fine gain control. At btock 991 , the present 
inventton performs the processing in an iterative manner ever a number of times (a programmabte nurrtier) to ref ne the 

30 calibration and adjustment When done, processing flews to btock 993 where the newly adjusted gain values for each 
channel (central and obscured) are stored in memory 1 102 (or otfier storage unit). 

It is appreciated that processes 970 and 980 can be performed by the camera system at any time. By performing 
the process 980 before imaging sessions of the gamma camera system, the preampfif ication gains of the PMT array 
are calibrated automatically such that the responses off the PMTs remain essentially the same as when they were cal- 

35 ibrated at the factory This will cause the scintillation detector to be caGbrated to match ttie response it had when ttie 
energy, linearity, and uniformity correction values were created. This match between the energy, linearity, and uniformity 
conection factors and the detector response improves image quality. 

In an alternative emtxxJiment of the present invention, in lieu of using a peak energy determination of the distribu- 
tions of processes 959. 969, 973. 985, and 987, an average or mean computation can be utilized. A "representative" 

40 data value of an energy or signal distrixition can therefore be ttie peak, average, mean, or similar function, of such dis- 
trbuti(^. 

In anotiier alternative emtxxliment of the present invention, the integrated channel signals for each obscured PMT 
are not measured to perform the gain correction for these PMTa In lieu of the integrated channel signals being 
recorded, the total energy of the events occurring witiiin the strip regtons 910 are recorded and an representative 

45 energy value is computed for ttiis distrixjtion for each obscured PMT. The above is performed for botti btocks 967 and 
983 such that the total energy of each ganrima event is recorded and compared. 

It is appreciated ttiat the data related to tiie integrated channel signals that is needed by ttie above procedure is 
gained via ttie computer system 1112 from ttie DEP 300 from bus 397 (see Rgure 2D). The total energy for a given 
event is supplied over bus 322 of DEP 300. In each case where channel information is gathered for gamma events 

50 detected wittiin spedfto XY regtons, the present invention provides a coincidence circuit for oonparing the XY coordi- 
nate of the event as output from bus 327 and 332 and ttie drcuit compares this coordinate value against the known 
re^pon (ag., area 570 or strip region 910). If a match occurs, then the desired data is sampled either via bus 322 fbr the 
total energy signal or via bus 397 for ttie specific channel data. Alternatively, ttiis specif to channel data can also be sup- 
plied directfy from buffer 325 of ttie DEP 300. 

55 

VARIABLE DYNAMIC COMPRESSION 

The present invention DEP 300 circuits (of botti detectors 80 and 80*) provides a programmabtoi variable dynamic 
compresston table 355 (Rgure 2D) so ttiat ttie energy response of a PMT can be acfusted using a programmable digital 
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function stored in menfx>ry rather tiian using dedicated analog hardvvare, as in the prior art In effect the dynamic com- 
pression tat}le 355 of the present inv^on Is utilized in place of the well known analog breakpoint driver circuits of the 
prior art The breakpoint driver circuits of the prior art are composed of analog circuits (including diode networks) for 
adjusting the PMT output. However, opposite to the present invention, the prior art circuits are not programmat>le or 
5 readily modlfiat)le. 

The dynamic compression tatAe 355 may be composed of programmable memory or static memory (ROM or 
PROM). If programmable memory is utilized, it may be composed of a number of weO known memory types, such as 
RAMorEEPROM. 

The dynamic corTf)ression table 355 of the present invention perfornrB two tasks: (1 ) it reduces or eliminates noisy 
10 signals from PMTs far from the point of the gamma interaction; and (2) its shapes the rentaining signals in order to 
obtain a more uniform spatial response across the detector. Digital implementation of this procedure allows the use of 
relatively simple linear preamplifiers thus redudng the complexity and cost of the analog electronics, and improving the 
long term stability of the system. The Dynamic compressk)n table 355 is utilized to alter the signal output from each 
PMT channel to remove nonlinearity so that when the signal summed with other altered signals, the resultant summa- 
16 tk>n signal is more linear. Refer to Rgure 1 4A whk:h illustrates the nomialized energy response of a PMT based on the 
distance from the PMT center at whk^h a gamma interactbn occurs. The response 652 is the signal t)efore correction 
by the dynamic compression table and response 650 illustrates the corrected or compressed response as output from 
tfie compression iogic of the DEP 300. Due to the characteristics and physk» of the PMT. the response is overty ele- 
vated for events that ocojr near the center of the PMT. therefore the conrpression to^ 
20 ingly. 

Refer to Rgure 14B. The overly elevated response of the PMT tut>e near the tut)e^ center poses a problem 
because when summed with other neigtixmng PMTs responses, the summation signal is not linear. The accuracy of 
ttie centroKl computatbn depends, in part, on the linearity of the summation signal as evkienced from the centroki com- 
putations as prevtously described herein. Figure 1 4B is a pfot of PMT response over PMT position for 5 PMTs in a par- 

2S tksjlar single axis. Response 654 is the summatfon of each of the 5 PMT signala F^re 14B illustrates that the 
summation signal 654 (sum of the PMT responses 655^) can be slightly nonlinear, and this nonlinearity tends to 
reduce the spatial accuracy of the centroid computation. The dynamic compression tatoUe 355 arx) subtractor drcuit 391 
provkies an adjustment to the input PMT integrated channel signal ever bus 352 to alk)w a more linear summation sig- 
nal. In aflM, after oonripensalkxit>y the dynanncoonrpressk^ 

30 stops. 

Rgure 15A illustrates a graph of the input energy response (normalized to global energy) 664 of a PMT channel 
and the desired normalized altered response 662 as a result of the dynamic compression procedure of the presmit 
invention. As sfiown, for moderate energy levels the higher the energy the oKxe the attenuation. Since tfie dynamic 
compresston procedure operates on normalized values, the dynamic compressk>n istHe 355 receives, as an input the 

35 gtobal energy over bus 322. The dynanuc compresston table 355 outputs, over bus 392. the value of the offset or differ- 
ence between signal 664 and 662 (for a given input signal over bus 352) mult9)lied by the gtobal energy (QE). There- 
fore, the output of the energy table 355 of the present invention is not normalized. The sut>tractor 391 then sut)tracts 
tfiis output of tat}le 355 from the uncorrected PMT channel signal stpplied over bus 352. The result is a corrected or 
compressed integrated cfiannel signal over bus 397. Figure 158 illustrates the response of the conversion utilized by 

40 the present invention for normalized input signals less than 15 percent of the max signal value. As shown, the corrected 
energy curve 662 is slightly higher in value than the input or uncorrected curve 664. 

In a particular embodiment, the coripression table 355 is 6 X 8 X 14 (32k) in size and is addr^^ 
of the global energy signals (322) and the 8 MSBs of the calibrated channel signal (352). Each 2 byte cell contains a 
14 bit offset value which is subtracted from the 10 bit value of the channel signal 352. The dynamic compresston table 

45 355 is caknjiated to 14 bits of preciston: 8 significant bits plus 6 fractional bi^ 

weighted suns to take on a greater number of values, reducing the quantization effects of the integer position calcula- 
tion. The output of the offset value from table 355 is 1 6 bits and the input to the sut>tractor from bus 352 are the MSB 
1 2 bits. The output from the compresston table 355 of tiie present invention is at a higher resolution than the signal over 
line 352. 

so An exemplary and general compression procedure utilized by an embodiment of the present invention is illustrated 
as the response shown in Rgure ISA and Figure 15B. Since the dynamto compresston is a non-linear function, the PMT 
integrated channel signals are scaled by the glok)al energy prior to dynamto compression. The scaling allows the 
dynanvc compression function to be independent of energy, which is an improvement over the non-linear amplifiers of 
the prior art 

ss Within this exemplary procedure, the integrated channel signal supplied over line 352 (Ep^t) is scaled, normalized, 
by the gtobal energy (QE) relative to the nominal peak global energy ^e) by: 

Ein = Epirt*(ge/GE) 
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Tlie value ge is oonslant and an exemplary value is suppGed below. Ihe scaM Ej„, is passed to the dynamic 
compression (lolloiO procedure, D, as shown further belcwL 

5 

Finally, the compressed signal is scaled back, denormalized. to its original level: 

S = Eft,-(QE/ge) 

10 where S is the signal that is output over line 387 from the subtractor drcuiL The exemplary general roll off procedure of 
the present invention is shown below as: 

D{E,^ = [SHi*e-^^"'""^''MSLo*[e-«^'^*"^^^1D .bias*S^ 

IS Where: 

SLo = l^Thresh * SP - Ejn (for Ejn <Lx>Thresh * SP) 
SLo =0(forEin>=LoThresh*SP) 
SHI = E|n - HiTresh * SP (for E^ > HiTresh * SP) 
20 SHI =0 (for Ejr,<= HiTresh *SP) 

According to the general procecbjre, SP is the nominal peak calbrated PMT signal and the each of the other parameters 
is expressed as a fractk>n of SR LoThresh Is the starting point of the low-end rollotf and HiTresh is the starting point of 
thehl^h-endroOoff. RLo cs the degrees of low-end rolloff. RHiisthedegreesof high-end roloff. Bias is the DC bias to 
25 be subtracted from the rolled off functkm. Preferred results have been achieved %inth the foHowing parameters (however 
each is programmable and adjustable within the present invention): 



30 



ge 


»180 


SP 


o200 


LoTresh 


»0.08 


HiTresh 


= 0.25 


RLo 


= 2.5 


RHi 


sO.35 


Bias 


= 0.05 



40 

in the computation of the preferred embocfiment of the present invention, the oompresskxi table 355 contains only tfie 
offset value that is sut>tracted from the cfiannel signal, e.g., tiie difference between lines 662 and 664 of Rgure 158. 
Therefore, the actual value stored in the lookup table 355 is computed based on the difference between ttie atX3ve pro- 
cedure and tfie chanriei input signal, times the gfobal eneigy GE. This is sficw^ 

45 

Data output over bus 392 = [ E,„-E J * [GE/ge] 

The output over bus 392 is then subtracted from the uncorrected energy of the particular PMT channel. Therefore, the 
procedure befow illuslrates the corrected or compressed signal: 

so 

Data oulput over bus 387 sEp^- Data output over bus 392 

The dynarnc compression circuit 355 of the present invention ensures that the above subtraction does not yieki a neg- 
ative number. Any negative number output is zeroed. 
SB The above procedure an be reduced and given t>elow assuming a normalized input signal from the range 0 to 1 00: 

= [ E^,*e-«' »'™^~ + ^^.(e-P^^oo . o.05 

where: 
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Ei„ = Ep^/GE 

TH=Ein-25. if Ejn < 25; and TH= 0. if Ej^ <=25 
TL = 8 - Ejn. if Ejn < 8; and Tb=0. if > s 8 
0.05 - small baseline offset 

5 

The input to the compression table 355 is GE (over bus 322) and Ep^t (over bus 352) for a particular PMT channel. The 
output of the table 355 is shown below. The actual value stored in the lookup tak)le is computed based on the difference 
between the above procedure and the normalized input, times the global energy GE. This is shown below: 

10 Data output over bus 392 = [ E ,„-E fj * GE 

The output ever bus 392 is then subtracted from the uncorrected energy of the particular PMT channel. Therefore, the 
procedure below illustrates the corrected or compressed signal: 

IS Data oufxjt over bus 387 = Ep^- Data output over bus 392 

The dynamic compression circuit 355 of the present invention ensures that the above subtraction does not yield a neg- 
ative number. Any negative number output is zeroed. 

It is appreciated that given sufficient memory size within drcuit 355 of the DEP 300, the data can be configured 

20 such that the global energy value GE is input along with the channel data of bus 352 and the compression table 355 
can then output the corrected value of the PMT channel data directiy over bus 387. In this embodiment ttie subtractor 
380 is not utilized. Also, the data of the dynamic compression table 355 of the present invention can be further opti- 
mized by modifying the procedures illustrated abova These modifications or optimizations may involve empirical data 
that is particular to a camera system or opeiating environment 

2s As shown in Rgure 15C, anottier dynamic compression procedure is illustrated that is stored in the dynamic com- 
pression table 355 and used in a similar manner with similar mechanisms as shown and discussed atxsva 

Given tiie above procedures, the actual data tables (addressed by GE and Ep„^ may readily be determined by one 
of ordinary skill in the art and are exemplary only. Because tfie table 355 is programmat)le, the level of resolution 
desired by any one embodiment is variable depending on the available size of memory for drcuit 355. According to the 

30 present invention, the data of ttie dynamic compressk>n table 355 is generated by data processor 1112 according to a 
programmable set of parameters and the compression procedure shown above and this information is ttien downloaded 
into the memory 355. If RAM is utilized, ttien the actual correction data can be stored in a disk storage 1 104 or ROM 
1 1 03 and tiien down loaded into the memory 355 before the gamma camera ^stem is used. Therefore, ttie data proc- 
essa 1 1 12 nriay store a variety of differem procedures and/a correction data tables and a u^ 

35 lor generating tables of correction data and for toading ttiem into memory 355 Ibr use in ttie DEP 300 as system cali- 
bratioa In ttiis altemativa ttie dynamic compression table 355 is extrerneiy flexible and pro^ammabla Alternatively, 
when the memory 355 is EEPROM, the data is programmed into ttie memory 355 before use and can be repro- 
grammed by a technician. 

40 AUTOMATIC BASELINE COMPENSATION 

Wrttiin the present invention, a sligtit DC offset (t>aseline offset) is maintained wittiin each preamplif teation stage for 
each PMT channel via a digital input to a DAC which controls inputs 212 and 214 of Figure 2C. This is ttie case for botti 
detector pair 80 and 80*. Shown in Figure 1 6 is an exervplary PMT analog channel output voltage signal 680 as sam- 

45 pled by ADC 236 (of Figure 2C). Due to ttie integration processes perfomned for an event, negative voltages of signal 
680 are unwanted, become dipped, and lead to inaccurate integration results. Due to drcuit drift of the amplifiers of dr- 
cults 280(0) to 280(54) (e.g.. due to temperature, current flow, etc.) the basefine of ttie voltage signal 680 as sampled 
t)y the ADC 236 may vary and, if left uncorrected, may dip t>elow zero. The t)aseline shift of an individual PMT is typically 
a fraction of a channel signat butthe sum of the shifte over an entire detector can move the giddal energy peak by sev- 

50 eral channela Baseline shifts also result in image registration variationa Drift may be compensated for Isy ttie use of 
expensive amplifier drcuits, but since a separate set of amplifiers is required for each channel, such a solution is not 
economtcalty practical. The present invention provides an automatic procedure and drcuitry lor baseline compensation 
in order to adjust for ttiis variable drift by measuring and adjusting the baseline voltage in real-time. 

In effect, ttie present invention inserts (via inputs 21 2 and 21 4) a base line compensation voltage amount called Vm 

55 in order to adjust ttie analog dgnal per channd. This is applicable for detector pair ^ Each channd has an 

independent Vm. This is shown in Figure 16. Line 682 shows true zero voltage and the signal 680 is a^usted by an 
anxxint Vm. After digitization and integration, ttiis offset value Vm must be subtracted from ttie total, resulting in ttie net 
channel signal. Basdine subtraction is peribrmed in ttie calibration lookup table drcuit 315 of ttie DEP 300 for each 
channel, and is based on ttie particular Vm inserted for ttiat channd. 
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The amount of voltage ofliset applied at 21 2 and 21 4 provided by the present invention varies* as discussed below. 
Hovvever, there is an ideal number Vm, that is computed based on the level of expected noise witfvn the dynamic range 
of the ADC 236. Since this dynamic range is within 1024 units, in one embodiment, and since the expected noise per- 
centage is roughly 8 percent, the digital value of Vm ^ approximately 78 units wHhin the dynamic range of the ADC 236. 

5 Of course this value will vary depending on the specific implementation of haidware utilized. This value is initially used 
as an offeet to control inputs 212 and 214 (a separate input is used for each channel) to offset signal 680. The value of 
Vm is determined such that signal 680 does not vary below zera Vm is also adjusted such that a large portion of the 
dynamic range of the ADC 236 is not consumed by the offset voltage. 

The actual value of the voHage offset as measured by the ADC 236 (when no event is present) may drift In the 

10 present invention, the digital processor 1 1 12 measures the signal output value for a particular channel that is not sut>- 
stantially receiving energy from a gamma event This value should be near the ideal voltage Vm. The sarrpled voltage 
is compared against the ideal voltage and is modified up or down (in a real-time feedback arrangement) to closely 
match the ideal voHage Vm. In this manner, the present invention compensates for the drift associated with the elec- 
tronics of the preamplificatibn stage for each channel by modifying the value of the offset voKaga This sample and 

IS adjust procedure is performed separately for each cfiannel. 

The present invention procedure for monitoring the analog signal of each channel and for modifying the baseline (rf 
this signal based on the ideal baseline value of Vm is described below. Generally, the present invention provides two 
alternate methods for baseUne Compensation: (1) a method used during periods of low count rate; and (2) a method 
used during periods of high count rate. The cotriputer system 1112 contains a baseline offset matrix of data values in 

20 memory 1 102. one for each channel, that records the current baseline offset voltage applied at Inputs 212 and 214. The 
values within this baseline offs^ matrix are updated by the present invention to conpensate for voltage drift of the 
anplifler stages of the circuits 280(0) to 280(54). Although described with respect to one detector, it is appreciated that 
the below procedure operates with respect to both detectors 80 and 80*. 

Ref^ to Figure 17A which illustrates the channel baseline compensation procedure 703 utilized t)y the present 

25 invention and perfonned by digital processor 1112. Procedure 703 is peribnned both respect to both detectors 80 and 
80'. At block 71 0. the present invention initially detennines if the count rate detected by the camera system is high or 
low. High Count rate is determined as a Coixii rate that exceeds 75.000 counts per second, but this court rate is exem- 
plary and is programmable within the present invention. Different procedures will be implemented tiy the present inven- 
tion depending on the camera'6 count rata At 715, the digital processor determines if the count rate is above the 

50 threshokj (e g.. 75.000 counts/see) and. if tme, block 760 is perfonned for event triggered baseline correction. If the 
overall count rate detected by the gamma detector is kwver than or equal to the threshokJ amount, then at block 720, 
the processor 1112 performs software triggered baseline confection. At the corrpletion of a correction procedure, the 
processing is exited via 799. It is appreciated that the f tow 703 of the present invention is repeated throughout the oper- 
atfonal duty cyde of the gamrna cainera induding periods in which the cainera is Mle a^ 

35 the camera is activeiy engaged in an imaging sessfon. 

Software Triggered. Refer to Rgure 17B which aiustrates the software triggered baseline correction procechire that 
is peribrmed during periods of low count rate. Rgure 17B illustrates the steps peribrmed by the present inventfon by 
block 720. At block 724, for each of the channels of the d^ector 80. the present invention inserts a software generated 
trigger that emulates an actual event trigger that would be d^ected over One 1 30 cf drcurt 1 00 of Figure 2A. For each 

40 channel (of both detectors), the digital processor 11 12 in effect triggers an integrator (e.g.. 238 or 240) to begin inte- 
grating over a pulse period. The integratfon is performed simultaneously for each channel at btock 726 (or altematively 
in stages) arvJ each channel should, kteally. integrate over the baseline voltage since statistically no event is expected 
to coincide with the software triggered pulse period. The software trigger is also called a false trigger since it is not ini- 
tiated t)ased on detection of a true gamma interaction. Since the software triggered Integratfon does not usually coin- 

45 dde with an event for the rnost part the integration peribrmed t>y the chm 

baseline voltage should, ideally, be very dose to Vm. This is the case because the channel shouki not be generating 
an output signal due to the expected (e.g.. statisticalty) absence of an event 

At block 728. the data generated from the integration of bfock 726 is recovered from tiie raw FIFO memory 310 (of 
DEP 300) by computer system 1 1 12. This memory contains the integrated result of each channel without subtraction 

so of tiie Vm offset (which is peribrmed t>y tiie calibration table 315). Optionally, at block 730. ttie data for each channel is 
divided by the pulse period to normalize the integration in tima This step is optional because If the time period of each 
integration is known, tiien tiiis time period multiplied by ttie ideal Vm value can be used as a reference by which ttie 
base line offoet value can be compensated. Either implementation is wrttiin the present invention. At block 732. the data 
for each channel is stored into a separate histogram that is maintained by computer system 1 1 1 2 in memory 1 1 02 for 

ss each channel. The histogram records the sampled baseline voltage (for a given channel) over tiie number of samples 
taken. At bfock 734, Vie above process is repeated for a next sample until a predetermined and programmable number 
of samples is taken. For each sample, tiie data associated witti each h^ograms are increased. An exemplary number 
of samples used by the present invention is 250 samples per channel. Therefore, a particular histogram constructed for 
each channel by btock 732 is composed of roughly 250 data pdntSw These histograms are maintained by the present 
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invention in RAM ntemory 1102 of the data processor 1112. Assuming 55 channels per detector, there needs Id be 
13.750 false counts per caIit)ration process of the present invention for process flow 720. 

At k)lock 736 of Figure 1 7B. the present invention utilizes the histogram created for a particular channel of a given 
detector to determine the average value of the sampled baseline offset amount This deterrranation may be made using 

5 a number of well known averaging, mean or weighted averaging procedures. Any nurTi>er of cfifferent procedures may 
be utilized consistent within the scope of the present Invention at block 736. At block 738. for the particular channel, the 
average Is compared against Vm (which is stored in memory by processor 1 11 2 and is constant for each channel). If 
the average for the particular channel is greater than Vm. then at bk)ck 742 the baseline offset of the particular channel 
which is stored in a baseline offset matrix within memory 1 1 02 (as applied via signals 21 2 and 21 4) is decreased by 

10 one unit and the matrix is updated and block 744 is entered. If the average for the particular channel is less than Vm. 
then at t3k>ck 740. the baseline offset of the particular channel (as applied via signals 212 and 21 4) is increased by one 
unit and the baseline offset matrix is updated and bkx:k 744 is entered. Otherwise, if the average value of the baseline 
is equal to Vm, then no baseline correction is required and tkxk 744 is entered drectly It is appreciated that there is a 
separate computer controlled pair of Inputs, 212 and 214. for each preampffication circuit 280(0) to 280(54). 

IS At bk)ck 744. the above procedure (to block 736) is repeated for each channel, so that the baseline corrections for 
each channel are adjusted (in a feecfl^ck arrangement) separately At block 746 the procedure is eodted and the histo- 
grams maintained in memory 1 1 02 for the above compensation procedure are reset TTierefbre. as can be seen above, 
the present inventton measures the channel voltage for each channel during periods of no expected events and records 
and adjusts the basefine voltage for the channels Due to drift, these voHages 

20 value Vm to prevent clipping (due to negative voltage). The present invention provides a measurement and compensa- 
tion procedure to effectively maintain the baseline voHages for each channel to the ideal value. Vm. This is performed 
using a feedb^ arrangement as shown in Figure 1 7B thus avoiding the use of expensive circuitry witNn the preampli- 
f k»tion stages 280(0) to 280(54). 

It is appreciated that the present invention, for bkxk 720, does not determine in advance if there is a real event 

25 occunring in coincidence with the software trigger generated at block 724. Statistically, it is rare (but possible) to have 
an event occur in coincidence with (and therefore interf^e with) the software triggered integration period during periods 
of low count rate. Should an event occur during bk>ck 726 of ttie present invention, that data point will be much larger 
than the expected value (due to the PMT responses) and will be effectively excluded statisticaUy from the average (or 
median) conputed Idr the channel at bkKk 736. In an alternative embodiment, tiie present invention may peribnn an 

30 additional step of excluding all data points from the histograms generated in block 732 that exceed a predetermined 
ttueshold in order to exclude data points cormpted by an event However, spven 250 samples (for instance), the number 
of data points corrupted by an event is very k3w and is negligible. 

Scintillation Triooered. Refer to Figure 17C and Figxre 17D which illustrate the procedure 760 utilized by the 
present invention to provide baseline compensation during periods of high count rates. Procedure 760 utilizes actual 

35 events to trigger the measurement processes and therefore is called event or scintillation driven baseline connection. 
The process begins at the detection of a gamma event at 764 and flows to 766. An 

over line 1 30. For each channel, the trigger pulse causes each integrator to integrate over the pulse period as described 
previously herein. Also at tkxk 766. the data from ttie integration is fbmvarded to ttie DEP 300. At block 768. the DEP 
300 detemrines which channel was ttie peak channel (e.g., ttie peak PMT) for ttie event by using circuit 320 (Rgure 2D) 
40 described at)ova 

Based on ttie peak PMT value (which is reported back to ttie digital processor 1 1 1 2). ttie present invention utilizes 
a memory circuit 800, as shown in Figure 18. to determine which PMTs are spatially near to ttie peek PMT and which 
PMTs are spatially far from ttie peak PMT. Fax PMTs are determined to be ttiose PMTs ttiat receive substantially no (or 
negligibly litfle) energy from ttie event and near PMTs are ttiose PMTs ttiat receive some energy (ag. , less ttian 0. 17%. 

45 etc.) from ttie event. As shown in Figure 1 8, ttiere is an entry for each PMT in column 805 and column 81 0 lists ttiose 
PMTs ttiat are near and colurnn 820 contains a list of ttiose PMTs ttiat are far. For each peak P 
far PMTs ttiat should receive subslantially no energy from tiie event For any given event detected at block 764 it is ttie 
far PMTs ttiat are sanpled by ttie present invention for baseline correction. The set of far PMTs for a given peak PMT 
address may be determined experimentally by measuring PMT responses for a kncwn peak PMT address wittiin a 

50 given PMT array. 

At block 770 of Rgure 1 7C, the present Invention addresses circuit 800 wrtti ttie peak PMT address (number) to 
determine ttie set of far PMTs as reported by column 820. At block 772. for ttie given event ttie present invention drects 
computer system 1 112 to recover ttie data associated witti each of ttie far channels as stored in ttie RAW FIFO memory 
310 of DEP 300. These far channels should report essentially an integration of ttie baseline voltage because ttiey 
55 received an insubstantial amount of energy from ttie detected event. Optionally, at bk)ck 774 ttie integration data of ttie 
droit 310 is divMed by ttie pulse period to namalize the data in time for each far channel. At tkxk 776, ttie present 
invention adds ttie value of ttie sanpled voltage to ttie histogram for each far PMT channel sampled. At block 778. ttie 
above procedure is repeated for the next sanf)le. This procedure continues until each channel obtains at least a prede- 
temnined and programmable number of data points in its histogram. An exemplary number of data points is 2S0 sanv 
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pies per channel. 

It \s appreciated that for each gamma event at 764, a different number and group of far channels will be sampled. 
Therefore, the number of histograms updated for each went will vary. The present invention therefore maintains a 
record of the channel having the least rumber of data points in its histogram and when fhas channel reaches the prede- 

5 termined number (ag.. 250). block 778 will exit to "A" as shown. 

Refer to Rgure 1 7D illustrating the remainder of process 760 of the present invention. From "A," the flow continues 
to block 780 where the present invention determines tiie average (ag., weighted average, median, meaa etc.) t>asefine 
voltage for each histo^am for all the channels. This process is anak)gous to the process of tkxM. 736 of the software 
triggered procedure 720 except each histogram of block 780 may contain a different number of data points depending 

10 on the spatial distnbution of the detected gamma events. At thck 782. the present invention determines for a particular 
channel if its average baseline voltage amount is larger than Vm. If so. then at ttock 786 the offset voltage stored in the 
baseline voltage matrix (and as asserted by signals 212 and 214) is deaeased for ttiat channel and block 788 is 
entered. If the average voltage is less than Vm for a partkailar channel, then at block 784 the baseline voltage as stored 
in the baseline voltage matrix is increased for that channel and \okxk 788 is entered. If the average baselnie offM volt- 

f 5 age fa the channel is equal to Vm, then no corection is required and block 788 is entered drectiy. 

At block 788. the present invention repeats from bkxk 780 for the next cfiannel until all baseline voHages of chan- 
nels of the detector head are compensated. At tkxk 790. the process 760 exists. It is appreciated that in the event that 
two gamma events are detected within a small window of time, one event may trigger the calibration sample (e.g., tkKk 
764) and the other event may interfere with the integration of the far PMTs and therefore increase or interfere with the 

20 sample data. In these statistically rare cases, as with the software triggered procedure, these corrupted data points are 
effectively excluded (or minimized) statistically via the averaging functions of the present invention. 

Therefore, as shown above, this embodiment of the present invention provides a mechanism and process that 
maintains the baseline offset voltage, per channel, to the kieal value Vm. Drift associated with the transistora of the 
preampliffoation stage of each channel is detected and corrected by using a feedback arrangement. By using the 

25 approach of the present invention, the preamplif ication stages of the detector circuitry m^ be implemented using less 
expensive (ag., less drift tolerant) circuitry 

It is appreciated that during periods of high count rate, the software triggered sampling procedure 720 of the 
present invention becomes less accurate as more events occur during the calibration samples. Therefore, the event 
driven procedure 760 is more accurate for high count rates. At kyw count rates, the event driven calibration procedure 

30 760 is not operational because the count rate becomes too fow, e.g., background radiation on the order of 40 counts 
per second is too low to perform event driven calitxration. Procedure 760 requires that each channel t>e exposed to at 
least 250 counts per baseline adjustment Therefore, the use of both software driven 720 and event driven 760 baseline 
compensation of the present invention is advantageous for high and low count environments. 

SS SWrrCHABLE PET AND SPECT M ODES OF OPERATION 

Figure 19 illustrates a process 1400 utilized t}y the present invention when switching between SPECT and PET 
imaging modes of operation. When operating within PET imaging modes, the patient is injected with a Ffouro Deoxi Gu- 
cose (FDG) or similar radiopharmaceutical prior to imaging emitting two 511 keV gamma rays. In SPECT moda the 

40 patient is injected with a TI-201 radtopharmaceutical prior to inriaging emitting a 140 keVgamnia ray. Process 1400 of 
the camera system of the present Invention starts at block 1405 wherein the system checks the user ff^^ 
(ag.. 1 106 of Rgure 5) or memory 1 102 to determine if a PET or SPECT imaging nrxxJe of operation is requested. If 
PET nrxxJe of operation is requested, then at block 1410. the camera system verifies to determine if the collimators are 
removed from the detectors 80 and 80*. Processing is suspended until the collimators are removed and an indication is 

45 displayed on unit 1 105 or 1065 to indk^ate the presence of the collimators. If the collimatora are removed, then the proc- 
ess contmues to btock 1415 wherein the camera system progranrs the trigger mode (ag.. line 1252 of the CTC 1050 off 
Rgure 2A3) to be in P ET imaging moda This will require coincidence trigger detection to generate a valid event trigger 
over lines 1240 and 1242. 

At block 1420, the computer system loads the PET integration (accumulation) interval (ag., approximately 320 ns) 
50 over line 235 of Figure 2C into register 1330 for both detectors 80 and 80*. At bfock 1425, the present invention directs 
the computer system to set or download a PMT cluster table having a smaller PMT cluster size (ag., 7 PMTS) into the 
PMT address table 335 of Rgure 2D for each indexed peak PMT This is accomplished by programming the PMT 
address table (see Rgure 8) so that the cluster sizes are smaller for each peak PMT. At bfock 1430. the present inven- 
tion adjusts the configuration of tiie detectors 80 and 80* such tiiat they are at a 180 degree configuration for imaging. 
55 At block 1435, the present invention perfonns a PET imaging sessfon across a number of projection angles with the 
detectors 80 and 80' in their 180 degree configuration. For each projection angle, within block 1435, data is collected 
by the present invention on an event by event basis according to the f tow diagrams shown on Rgure 4A and Rgure 4B. 
Within PET moda coincidence detection is perfonned between detectors 80 and 80* for the events. Collected data is 
recorded into the aoqiosition computer system 1055 (Figure 1). It is qspreciated that in PET mode; the energy value, Z, 



27 



EP0747 730A2 



output by the detectors in response to an event is the sum of the energy detected by the cluster used 1o locafize the 
event and not the sum of all the PMTS in the detector as is sometimes used in SPECT imaging. The collected data is 
then imaged by inr^age processa 1060 and PET reconstruction Is performed on the projection data at block 1440. PET 
reconstruction procedures on projection data are well known and are not discussed in detail herein. It is appreciated 
5 that at btock 1 440. attenuatk>n connection maps can be applied to the emisston data to conrect for nonuniform attenua- 
tk>n. 

At block 1485, a user-selected image from the PET reconstruction can be displayed on display screen 1065 or cfis- 
piay 1 105. Since PET imaging data is collected without a collimator, the resolution of PET images is generally of higher 
(^laOty over SPECT mnages. It is appreciated that the processing of blocks 1410 - 1430 can occur in any order and the 

10 order presented is exemplary only. 

With respect to Rgure 19. if SPECT mode of operation is requested at block 1405. then at block 1450. the camera 
system verifies to determine if the collimators are present on the detectors 80 and 80*. Processing is suspended until 
the collimators are installed and an indication is displayed on unit 1 105 or 1065 to indicate the absence of the collima- 
tors. If the oolBmators are installed, then the process continues to tkxk 1 455 wherein the camera system programs the 

15 trigger mode (e.g.. line 1 252 of the CTC 1050 of Rgure 2A3) to be in SPECT ima^ng mode which does not require 
coincidence detection for valid event trigger generation. 

At block 1460. the computer system loads the SPECT integration (accumulation) interval (ag.. approximately 840 
ns) over line 235 of Rgure 2C into register 1330 for both d^ectors 80 and 80*. At block 1465. the present invention 
directs the computer system to set a relatively larger PMT cluster size (ag.. larger than 7, for instance, 19 PMTS) into 

so the PMT address table 335 of Rgure 2D for each Indexed peak PMT. This is accomplished by programming the PMT 
address table (see Rgure 8) so that the cluster sizes are relatively larger for each peak PMT At block 1470. the present 
invention adjusts the configuration of the detectors 80 and 80' such that th^ are at a 1 80 degree configuration for imag- 
ing. At bkKk 1 475. the present invention performs a SPECT imaging session across a number of projection angles with 
the detectors 80 and 80' in their 180 degree configuratton. For each projection angle, within bkxdc 1475. data is collected 

25 by the present inventkxi on an event by event basis according to the ftow cfiagrans shown on Figure 4A and Figure 4B. 
Within SPECT mode, coincidence detection Is not performed. Collected data is recorded into the acquisition computer 
system 1055 (Figure 1). The collected data is then imaged by image processor 1060 and SPECT reconstruction is per- 
formed on the projection data at bfock 1480. SPECT reconstructkMi procedures on projection data are well known and 
are not discussed in detail herein, ft is appreciated that at bfock 1480. attenuation correctfon maps can be applied to 

30 the emissfon data to con-ect for nonuniform attenuation. 

At bfock 1485 a user-selected image from the SPECT reconstructfon can be displayed on display screen 1065 or 
display 1 105. It is appreciated that the processing of bfocks 1450 - 1470 can occur In any order and the order presented 
is exemplary only. Using the procedure 1400. the present invention offers a nuclear camera system automatically 
switchable and optimized between SPECT and PET imaging capabilities. 

35 The preferred embodiment of the present inventfon. a dual head nuclear camera system automatically switchable 
between SPECT and PET modes of operation is thus described. While the present inventfon has been described in par- 
ticular embodiments, it should be appreciated that the pres^ invention should not be construed as limited by such 
embodiments, but rather construed according to the befow claims. 

40 Claims 

1. A nuclear camera system comprising: 

a first detector and a second detector for detecting sdntmatton events, sakJ fvst detector and saki second 

45 detector coupled to an acquisition conputer system; 

a input mechanism for selecting between a PET imaging mode and a SPECT imaging mode; 
a switchable trigger detection unit that is automatically switchable between coincidence detectfon when said 
PET imaging mode is selected and non-coinckJence detection when sakf SPECT imaging mode is selected; 
said switchable trigger detectfon unit covpH&i to supply a first valid event trigger signal to said first detector 

so causing saki first detector to record a sdntillatfon event and coupled to supply a second valid event trigger sig- 

nal to said second detector causing said second detector to record a sdntiDatfon event 

2. A nuclear camera system as described in Claim 1 wherein: 

said first detector and saki second detector incfivfoually comprise event detection units for transmitting incfi- 
55 catfons of scintillation events to sakJ smtchable trigger detectfon unit; and wherein 

saki 8witcfiat)le trigger detection unit generates sakJ first and second valid event trigger signals upon detec- 
tion of two indications of sdntillation everts occurring within a coinckJence window, one in saki first detector and 
one saki second detector, provkied saki PET imaging mode is selected. 
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3. A nudear camera system as described in Claim 2 wherein each event detection unit of said event detection units 
correspond to a zonal area of said first or second detector and each event detection unit generates an indication of 
sdntiUation events occurring within a corresponding zonal area of said first or second detector. 

5 4. A nudear camera system as descrbed In Claim 3 wherein said event detection units utilize constant fraction dis- 
criminators to generate said incfications of sdntillation events. 

5. A nudear camera system as described in Claim 1 wherein: 

said first detector and said second detector incfividually comprise event detection units for transmHting inJ- 
10 cations of sdntillation events to said switchable trigger detection unit: and wherein 

said switchable trigger detection unit generates said first valid event trigg^ signal upon receiving an indica- 
tion of a sdntillation event occurring within said first detector and generates said second valid event trigger signal 
upon receiving an indication of a sdntillation event occurring within said second detector, provided said SPECT 
imaging mode is selected. 

IS 

6. A nudear camera system as descritDed in Claim 5 wherein each event detection unit of said event detection units 
correspond to a zonal area of said first or second detector arxJ each event detection unit generates an indication of 
sdntillation events occurring within a corresponding zonal area of said first or second detector. 

20 7. A nudear camera system as described in Claim 6 wherein said event detection units utilize leading edge cfiscrimi- 
nators to generate said indications of sdntillation events. 

& A nudear camera ^stem as described in Claim 1 wherein: 

said first detector comprises PET event detection units for transmitting indications of sdntillation events 
25 occurring within said f irst detector to said switchable trigger detection unit: 

said second detector comprises PET event detection units for transmitting indications of sdntillation events 
occurring within said second detector to said switchable trigger d^ection unit; and wherein 

said switchable trigger detection unit comprises a coincidence detector coupled to said indications of sdn- 
tillation events occurring within said first detector and coupled to said indications of sdntfllation events occurring 
30 within said second detector, said coincidence detector for generating a cdnddence signal upon cdnddence of said 
sdntiDation events occurring within said first detector and said incfications of scintillation events occurring within 
said second detector. 

9. A nudear camera system as desabed in Claim 8 wherein: said first detector comprises SPECT event detection 
35 units for transmitting indications of scintillation events occurring within said first detector to said switchable trigger 
detection unit; and 

said second detector comprises SPECT event detection units for transmitting indications of sdntillation 
events occurring within said second detector to said switchable trigger detection unit. 

40 1 0. A nudear camera system as described in Claim 9 wherein said switchable trigger detection unit further comprises: 

afirst multiplexer receiving a select signal and receiving an input from said SPECT event detection units of said 
first detector and receiving an input from said coincidence detector; 

a second multiplexer receiver select signal and receiving an input from said SPECT event detection units of 
45 said second detector and receiving an input from said coincidence detector; and 

wherein said first multiplexer outputs said first vafid event trigger signal to said first detector and wherein said 
second nwKiplexer outputs said second valid evem trigger signal to 8^ 

11. For a channel drcuit associated with a photomultiplier of a sdntDlation detector, an accumulator drcuit comprising: 

50 

a programmable register for containing a programmable integration period; 

an analog to digital converter for receiving a first signal representative of energy detected by said photomulti- 
plier as a result of a sdntillation event and for digitizing said first signal to generate a digital signal therefrom; 
a first integrator for integrating said digital signal over said programmable integration period to arrive at a result; 
55 a memory drcuit for containing said result of said first integrator. 

12. An accumulator drcuit as described in Claim 1 1 further comprising a programnrting unit for programming a first inte- 
gration period into said programmable register such that said f irst integrator is optimized for periods of higher count 
rale and for progamining a second Integration period bito said programmable register 
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is optiinzed for periods of lower count rate. 

13. An accumulator drcurt as described in Claim 1 1 furttier comprising: 

a progfEunming urat for programming said programmable integration period such that said first integrator is opti- 
mized for a PET operational mode. 

14. An accumulalor drcuit as described in Qaim 1 3 further comprising: a programming mit for programming said pro- 
grammable integration period such that said first integrator is optimized for a SPECT operational nvxie. 

15. An accumulator drcurt as described in Claim 12 wherein a received trigger signal initiates an integration interval 
and further comprising: 

a counter drcurt for counting an intenml from said trigger si^ial; 

a comparator drcuit cotpled to said counter drcuit and coipled to said programmable register, said compara- 
tor drcuit generating an output signal when said counter drcuit reaches said programmable integration period; 
and 

said memory drcuit for latching in said result upon generation of said output signal from said comparator dr- 
cuit 

16. In a nudear medicine camera, an apparatus for generating a photomultiplier duster, said apparatus comprising: 

at least a pair of Scintillation detectors for reporting a pair of events occurring in cdnddence for PET imaging 
and for reporting events individuany for SPECT imaging, each of said pair of sdntillation detectors conrprising: 

an array of photomultipliers wherein individual photonrultipliers generate channel signals, 

integration drcuitry integrating said channel signals from each sdntillation detector responsive to a 

gamma event and generating integration resuHs therefrom; 

peak drcuitry for determining a peak photomultiplier based on said integration resuHs and for generating 
a signal indicative of said peak photomtMpiier. and 

duster d rcuitry addressed by sakJ peak drcuitry and responsive to said signal indicative of said peak pho- 
tonrtittiplier for generating a photontiltiplier duster associated with said peak photomultiplier, saki duster 
drcuitry containing separate pfiotonrultiplier dusters for indivkiual photomultipliers of said array of photo- 
multipliers and wherein saM cluster drcuitry also generates a cluster type signal responsive to said signal 
indicative of saki peak photonrultiplier wherdn said duster type sig^ 
photomultiplier duster. 

1 7. An apparatus as described in Claim 1 6 wherdn saki photomultiplier duster comprises a set of photonrultiplier indi- 
cators and further comprising: 

a sequencer circuit coupled to saki duster drcuitry. saki sequencer drcuit for addressing saki duster drcuitry 
such that indivkiual photomultiplier indk»lors of said pholornultipner duster are sequentially output from saki 
duster drcuitry; and 

weight drcuitry responsive to (I) saki photomultiplier indteators and (2) saki duster type signal, saki weight cir- 
cuitry for generating indivkiual coordinate weight value signals associated with individual photonuttipHers of 
saki photomultiplier cluster. 

18. An apparatus as desaibed in Claim 16 wherein said duster type signal indicates between a normal duster type, 
an edge duster type and a comer duster type of saki photomultiplier cluster. 

19. An apparatus as described in Claim 16 wherein saki photomultiplier duster comprises a set of photomultiplier indi- 
cators and further comprising a sequencer drcuit coupler to saki duster drcuitry. saki sec^encer drcuit for 
addressing saki duster drcuitry such that individual photomultiplier indicators of saki photomultiplier duster are 
sequentially output from said duster drcuitry. 

20. An apparatus as described in Claim 19 wherein saki cluster drcuitry comprises a programmable memory drcuit 
addressed by saki signal indicative of saki peak photomutt^plier and saki sequencer drcuiL 

21 . An apparatus as described in Clakn 1 7 wherein saki photomultiplier duster comprises a set of photomultipGer incfi- 
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catOTS and further conrprising: 

a buffer storage for storing integrated results of said photomultiplier duster, said tx/ffer storage addressed bf 
photomultiplier indicators output from said duster drcuitry: and 
5 c^itrotd computation drcuitry coupled to receive output integration results from said txjffer storage and cou- 

pled to receive said individual ooorcfinate weight value signals assodated with individual photomuKpliers of 
said photomultiplier duster, and in response thereta said centroid conputation drcuitry for generating a two 
dimensional spatial cooidinate of said gamma event. 

10 22. An apparatus as described in Claim 16 further comprising a user selectable operational mode t}etween said PET 
imaging and said SPECT imaging wherein said photonrultplier duster that is output from said duster drcuitry is 
responsive to said operational mode. 

23. An apparatus as described In Claim 22 wherein said duster drcuitry outputs photomultiplier dusters in sets of 
IS seven or less photomuHipliefs each, provided said operational mode is set for said PET imaging and wherein said 

cluster drcuitry oulputs photomultiplier dusters in sets of more than seven photomultipliers each, provided said 
operational mode is set for said SPECT ima^'ng. 

24. In a nudear medidne camera, an apparatus for optimizing said camera for either PET or SPECT ima^ng modes, 
20 said apparatus comprising: 

a switchable operational nrKXIe between said PET imaging and said SPECT imaging nxxies; 
at least a pair of scintillation detectors for reporting a pair of events occuning in coincidence for said PET Imag- 
ing mode and for reporting events individually for said SPECT imaging mode, said pair of sdntillation detectors 
25 for rotating akxxit a patient and for collecting projection image data under said PET or SPECT imaging modes: 

programmable integration drcuHry coupled to channel signals output from photonuHiplier tubes of said pair of 
sdntillation detectors, said programmable drcuitry for integrating said diannel signals over a first integnrtion 
period or over a second integation period in response to said operational mode; 

a programmatiie trigger ctrcuit coupled to event detection logic within said pair of sdntillation detectors, said 
30 programmable trigger drcuit for generating valid event trigger signals based on coincidence between said pair 

of sdntillation detectors or not bdsecH on coincidence k)etween said pair of sdntillation detectors in response to 
said operational mode; arxJ 

a programmable photomultiplier duster drcuit responsive to values over saM channel signals, said program- 
mat)le photomultipGer cluster for utilizing a duster of a f ffst size fbr localizing detected events and fbr utifizing a 
3s duster of a second size for localiznig detected events in response to said operational moda 

25. An apparatus as described in Claim 24 wherein fast integration period is utilized for said PET imaging mode and 
wherein said second integration period is utilized fbr said SPECT imaging mode arxJ wherein said first integration 
period is shorter in length over said second integration period. 

40 

26. An apparatus as described in Claim 25 wherein said programmable trigger drcuit utilizes coincidence between said 
pair of sdntillation detectors for generating said valid event trigger signals when in said PET imaging mode and said 
programmatjie trigger drcuit does not utilize coincidence t>etween said pair of sdntillation detectors for generating 
said valid event trigger signals when in said SPECT imagbig nrxxJe. 

45 

27. An apparatus as descrit>ed in Claim 26 wherein said first size of said programmable photomultplier cluster is seven 
photomultipliers or less and is utilized when in said PET imaging mode and wherein said seoorxJ size is more tiian 
seven photomultipliers and is utilized when in said SPECT imaging mode. 

50 



55 



31 



EP0747730A2 




32 



EP0747730A2 




llOd 



1132 




FIG.2A-1 



33 



EP0747730A2 



Ao- 
Al7- 



A30— 



A24 — 



A43-»- 



A37— 



A54— 



EVENT DETECTION 



EVH^ DETECTION 



EVENT DETECTION 



HOC- 



EVENT DETECTION 



/ 

Detector go 



llOd' 



IQQ" 



M26b 




>128b 




All- 
A30— 



A24 — 



A43-*- 



llOa' 



EVENT DETECTION 



.l^Ob* 



EVENT DETECTION 



110c' 



A37— 
A54-^ 

/ 

Detector 30' 



EVENT DETECTION 



llOd" 



124b 



M26b 




-128b 



FIG.2A-2 



34 



EP0747730A2 



1050 



1130 



1130' 

> 



p->^1125 



1132" 




1135 



1132 



1310 



1240 



Valid Event Trigger 



1252 



-< PET/SPECTMode 



W1320 



-^i;- ^ Valid Event Trigger 

1242 



FIG.2A-3 



35 



EP0747730A2 



Ao 



R 



A54. 



R 



144 

FIG.2B 



36 



EP0747730A2 




EP0747730A2 




38 



EP0747730A2 



- I 



CM 

s 



CO 



2^ 



21 £i 

X K 



o 




CO 
LL 



39 



EP0747730A2 



460 



START 

31 



BASEUNE 
OFFSET PER 

CHANNEL 



GAIN 
ADJUSTMENT 
PER CHANNEL 



J. 



DIGITIZE 
SIGNAL 



462 

464 
< 



TRIGER 
DETECTION 



466 



468 



INTEGRATE ON TRIGGER 
FOR EACH CHANNEL AND 
STORE 



470 



1 



INTEGRATE ON TRIGGER 
FOR EACH CHANNEL AND 
STORE 



472 



474 



TRANSFER 
CHANNEL 
DATATODEP 



478 



462 



DATA 
PROCESSOR 
SAMPLE RAW 
CHANNEL DATA 



476 



REMOVE OFFSET 
ADJUST FINE 
GAIN PER 
CHANNEL 



STORE ALL 
CHANNEL DATA 
IN BUFFER 



460 



DETERMINE 
reAKPMT 
ADDRESS 



COMPUTE 
GLOBAL 
ENERGY 



484 



DETERMINE CENTROID 
TYPE,SHAPE,AND 
PMT ADDRESS FOR EACH 
PMT OF CENTROID 




486 



FIG.4A 



40 



EP0747730A2 




BESET 
CENTROID 
LOGIC 



4B6 



GET NEXT PMTT 
ADDRESS OF CENTROID 

490 



DETERMINE 
WxANDWy 



492 



DETERMINE 
ENERGY 
(DYNAMIC 
COMPRESSION) 
493 



NO 



MULTIPLY AND 
ACCUMULATE 
x.y, ENERGY 
495 




DIVIDE BY 
LOCAL ENERGY 

497 



OUTPUT (x.y) 
COORDINATE 



498 



c 



RETURN 



T 



RETREVE 
INTEGRATED 
SIGNAL DATA 
FROM BUFFER 
494 



FIG.4B 



41 



EP0747730A2 




42 



EP0747730A2 




-"1- 
z 



43 



EP0747730A2 




44 



EP0747730A2 




45 



EP0747730A2 



O 
u 



i4 

m 
iir 



c 










1 












• 
• 
• 












































CM 












































O 






















cn 






















CD 






to 






















CO 
r— 
















iD 






cn 




STOF 












in 






01 
r— 




N 
CO 




STOP 














o 

CM 




Ol 
CO 




CO 

yt 








tn 






00 




CO 
CM 




CM 








CM 










CM 
CM 




tn 
in 












m 
in 




tn 
tn 




UJ 






NORMAL 




EDGE 




CORNER 










o 


* * • 


1^ 


• • • 


00 
CO 


• • • 


iD 


• • • 


•<*■ 

tn 





CO 

O 

Li. 



5" ^ 
DL m 



5 



46 



EP0747730A2 




47 



EP0747730A2 




48 



EP0747730A2 




49 



EP0747730A2 



570 




140 kev GAMMA EVENT 

FIG. 12 



50 



EP0747730A2 




WTTH COUMATOR REMOVED, MEASURE 
THE TOTAL ENERGY RESPONSE FOR 
GAMMA INTERACTIONS DETECTED 
WIDER REGIONS 570 ASSOCIATED 
WTTH EACH PMT IN THE DETECTOR 



I 



COMPARE THE MEASURED TOTAL 
eiERGY PEAK VALUE WfTH NOMINAL 

TOTAL ENERGY PEAK VALUE FOR 
EACH PMT AND DVLCULATE NEW PMT 
GAM VALUES ^ 



ADJUST GAIN VALUES OF THE 
CHANNELS BY CHANGING GAM 
VALUES IN THE PREAMPUFiCATION 

CnCUITS 280(1) AND OR THE 
CALIBRATION TABLE 315 FOR EACH 
PMTCHANNa 961 



NEXT 




STORE NEW GAIN VALUES IN MEMORY 
OR SIMILAR STORAGE DEVICE 

965 



I 



WITH COUMATOR REMOVED, MEASURE 
THE PMT RESPONSE FOR GAMMA 
INTERACTION THAT OCCUR WITHIN 
STRIP AREA 910 FOR ALL OBSCURED 
PMTS ^ 



I 



STORE SPECTRUM PEAK VALUE FOR 
EACH OBSCWED PMT IN MEMORY OR 
SIMILAR STORAGE DEVICE 




FIG.13A 



51 



EP0747730A2 






MEASURE THE TOTAL ENERGY 
RESPONSE FOR GAMMA 
INTERACTIONS DETECTED UNDER 
REGIONS 570 AS50CIA I cU Wl 1 n tAUM 
PMT IN THE DETECTOR 


1 ► 








COMPARE THE MEASURED TOTAL 

tiMcnuiT rCniX vnLUC Vflin INUmllMfM. 

TOTAL ENERGY PEAK VALUE FOR 
EACH PMT AND CALCULATE NEW PMT 
GAIN VALUES 




1 

T 




ADJUST GAIN VALUES OF THE 
CHANNELS BY CHANGING GAIN 
VALUES IN THE PREAMPLIFICATION 

CIRCUITS 2B0(i) AND OR THE 
CALIBRATION TABLE 315 FOR EACH 
PMT CHANNEL 075 


NEXT 


^^nEWOE until\ 



977 



DONE 



STORE GAIN VALUES IN MEMORY FOR 
EACH PMT IN MEMORY OR SIMILAR 
STORAGE DEVICE 

979 



RETURN 



FIG.13B 



52 



EP0 747730A2 






WITH COlilMATOR INSTALLED. MEASURE 

THE TOTAL ENERGY RESPONSE FOR 
GAMMA INTERACTIONS DETECTED UNDER 
REGIONS 570 AS ASSOCIATED WITH EACH 
NONOBSCURED PMT M THE DETECTOR^, 

Sol 


1 ^ 










WITH COLUMATOR INSTALLB) MEASURE THE 
PMT RESPONSE FOR GAMMA INTERACTIONS 
THAT OCCUR WITHIN STRIP AREA 910 FOR 
ALL OBSCURED PMTs 

983 






r 




COMPARE THE MEASURED TOTAL ENERGY 
PEAK VALUE WITH NOMINAL TOTAL ENERGY 
PEAK VALUE FOR NON OBSCURED PMTs 
AND CALCULATE PMT GAIN VALI£S FOR 

THOSE PMTs ggg 










COMPARE THE SPECTRUM PEAK VALUE 
MEASURED IN BLOCK 983 AGAINST THE 
VALUE STORED IN BLOCK 969 FOR EACH 
OBSCURB) PMT AND CALCULATE NEW GAIN 

987 






r 




ADJUST PMT GAIN VALUES BY CHANGING 
GAIN VALUES M THE PREAMUFICATION 
CIRCUITS 280(n ANUORTHE CALIBRATION 
TABLE 315 FOR EACH PMT CHANNEL 

969 





^991^ 



DONE 



STORE NEW GAIN VALUES IN MEMORY OR 
SIMILAR STORAGE DEVICE 

993 



REIURN 



] 



FIG.13C 



53 



EP0747730A2 



100 -r 




DISTANCE 
FROM PMT CENTER 

FIG.14A 




-2-10 1 2 X 



FIG.14B 



54 



EP0747730A2 




55 



EP0747730A2 




0 20 40 60 80 100 



ENERGY NORMALIZED 

FIG.15C 



56 



EP0747730A2 




57 



EP0747730A2 



zoa 



PERFORM 
SOFTWARE DRIVEN 
(TRIGGERED) 

BASELINE 
CORRECTION 




DETERMINE 
SYSTEM COUNT 
RATE 

710 




YES 



PERFORM 
EVENT DRIVEN 
(TRIGGERED) 

BASELINE 
CORRECTION 



760 



RETURN 



799 



FIG.17A 



58 



EP0747730A2 



C 



ENTER 



722 



5 



FOR EACH CHANNEL. 
MSERT SOFTWARE TRIGGER 

724 



FOR EACH CHANNEL. 
INTEGRATE OVER PULSE PERIOD 

726 



RECOVER CHANNEL DATA FROM 
RAW FIFO FOR ALL CHANNELS 
728 



FOR EACH CHANNEL 
DIVIDE BY PULSE PERIOD 



730 



STORE DATA FOR EACH 
CHANNEL IN INDIVIDUAL 

HISTOGRAM 732 



NEXT SAMPLE 




END SAMPLE 



DETERMINE AVERAGE 

FOR CURRENT 
CHANNEL HISTOGRAM 
736 



NO 



INCREASE BASEUNE 
OFFSET 

740 




(Optional) 



FIG.17B 



YES 



DECREASE 
BASEUNE OFFSET 

742 



NEXT CHANNEL 




END CHANNE 



RETURN 



748. 



59 



EP0747730A2 



760 




EVENT 



s ^ 






FOR 
EACH CHANNEL 
INTEGRATE OVER 
PULSE PERIOD 766 






r 




DETERMINE PEAK PMT 
CHANNEL 

768 







DETERMINE SET OF 
FAR CHANNELS 

770 



RECOVER RAW FIFO 
DATA FOR FAR 
CHANNELS 
772 



FOR EACH FAR 
CHANNEL DIVIDE BY 
PULSE PERIOD 

774 



GENERATE HISTOGRAM 
FOR EACH FAR 
CHANNEL 

776 



NEXT SAMPLE 




FIG.17C 



EP0747730A2 




DETERMINE 
AVERAGE FOR 
CURRENT 
CHANNEL 
HISTOGRAI^gp 




INCREASE 
BASELINE OFFSET 






DECREASE 
BASELINE 
OFFSET 




784 








786 






^ 


^ 







NEXT CHANNEL 
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SET TRIGGER 
MODE = PET 



1415 
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SET SMALL PMT 
CLUSTER SIZE 
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SET LARGE PMT /"^ 
CLUSTER SIZE 
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ADJUST 
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DISPLAY 
SELECTED IMAGE 
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RETURN 



FIG. 19 



63 



